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ENGINEERING  CONSIDERATIONS  ON  THE  SELF-ENERGIZING 
MAGNETOPLASMADYNAMIC  (MPD)-TYPE  FUSION  PL  ASMA  THRUSTER 


L  INTRODUCTION 

During  this  initial  period  of  the  three-year  project,  various  fusion  propulsion  concepts,  and  in 
particular  a  plasma  focus  device  and  its  analysis,  have  been  studied.  Based  on  the  identified  key 
operating  parameters  ( e.  g.,  capacitor  voltage,  capacitance,  inductance,  electrode  dimensions, 
plasma  pinch  current  and  temperature,  etc. )  of  existing  plasma  focus  devices,  feasibility  studies 
of  the  dense  plasma  focus  (DPF)  device  as  a  fusion  propulsion  thruster  have  been  performed. 

The  initial  findings  indicate  that  the  dense  plasma  focus  device  provides  attractive  plasma 
conditions  leading  to  the  ignition  for  DPF  as  a  space  propulsion  thruster,  though  uncertainties 
remain  in  tbe  validity  of  scaling  laws  on  capacitor  mass  at  high  current  beyond  1 MA.  Both 
conventional  and  spin-polarized  D-3He  fuels  were  studied  for  parametric  analysis  of  DPF  as  a 
propulsion  system.  The  major  portion  of  this  report  consists  of  the  parametric  study. 

Ever  since  Robert  Goddard,  Theodor  von  Karman,  and  Werner  von  Braun  proposed  concepts 
of  rocketaiy,  various  rocket  engines  have  been  promoted  as  the  main  vehicle  for  space 
exploration.  On  July  20, 1989,  President  Bush  announced  the  new  Space  Exploration  Initiative 
to  ruum  to  the  Moon  and  then  to  put  a  manned  presence  on  Mars.  The  goal  of  this  current 
pr  ject  is  consistent  with  the  Earth-Mars  mission,  with  the  emphasis  placed  on  a  fusion-powered 
propulsion  device. 

Conventional  propulsion  devices  based  on  chemical  energy  have  been  used  successfully  for 
early  space  explorations  with  firing  duration  ranging  from  seconds  to  a  few  hours.  Chemical 
propulsion  has  been  used  for  the  Saturn  5  (the  same  rocket  that  was  used  successfully  for  the 
Moon  Mission  in  1969),  the  space  shuttle,  interplanetary  probes,  and  the  missile  defense  system. 
The  best  current  chemical  engines  can  provide  specific  impulse  values,  Ijp,  (i.e.,  the  ratio  of  the 
rocket  thrust  to  the  propellant  weight  flow  rate )  approaching  500  seconds  [1], 

The  Nuclear  fission  engines  utilize  both  solid-  and  fluid-core  (liquid,  gas,  and  plasma) 
reactors.  Graphite  solid-core  reactors  operating  at  2000  to  2700°K  have  been  well 
demonstrated,  though  the  fluid-core  reactors  are  still  being  investigated  very  actively.  Typical 
specific  impulse  for  the  nuclear  fission  engine  ranges  from  600  to  1 100  sec  [1-2].  Other  rocket 
engine  concepts  (e.g.,  ion,  solar  thermal,  etc.)  have  also  been  considered.  Performance  of  220- W 
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and  10-kW  class  xenon  ion  thrusters  have  been  investigated  recently  in  Europe  and  at  NASA 
laboratories,  respectively  [3].  Rare  gas  propellants  (xenon  or  argon)  rather  than  mercury  are 
being  used  in  ion  thrusters  to  avoid  any  spacecraft  surface  erosion  problems.  Solar  powered 
electric  propulsion  and  solar-pumped  lasers  for  space  applications  are  also  being  pursued  [4]. 
Typical  performance  parameters  for  different  rocket  engine  types  are  summarized  in  Table  1. 


Tabte  1.  Rang**  of  Typical  Performance  Parameter!  hr 
Several  Different  Rocket  Engine  Type*  l  Rah.  14  J 

Specific  Maximum  Thruat-to-  Specific  Typical 

Engine  type  Impulse  (sec)  Temperature  Weight  ratio  Duration  Power  (hpflb)  Working  Pluid 

CK) 


Chemical  (liquid) 

300  to  460 

4,500  to  7,800 

10'a  to  100 

Seconds  to 
a  few  hours 

0.1  to  1,000 

Hj  toOg 

Chemical  (solid) 

200  to  310 

4,500  to  7,500 

10-*  to  100 

Seconds  to 
minutes 

0.1  to  1,000 

fuel  and 
oxidizer 

Chemical  (hybrid) 

200  to  400 

4,000  to  7,500 

10"*  to  100 

Seconds  to 
minutes 

0.1  to  1,000 

fuel  and 
oxidizer 

Nuclear  fission 

600  to  1,100 

5,000 

10-*  to  30 

Seconds  to 
a  few  hours 

0.1  to  1,000 

Hj 

Radioactive 

Isotope  Decay 

400  to  700 

2.200  to  3,000 

10'*  to  10-J 

Daya 

0.001  to  0.01 

H, 

Arc  Heating 

400  to  2,000 

10,000 

10-4  to  lO-* 

Days 

0.001  to  1 

H, 

Ion 

5,000  to 
25,000 

to-5  to  ltr* 

Months 

0.001  to  1 

mercury, 
xenon,  argon 

Solar  thermal 

400  to  700 

2,500 

10'*  to  to-1 

Days 

0.001  to  1 

h2 

The  chemical  rockets,  and  to  some  extent  also  the  nuclear  fission  rockets,  have  relatively  low 
values  of  specific  impulse,  relatively  low  engine  weight,  a  very  high  thrust  capability,  and, 
therefore,  high  acceleration  and  high  specific  powers  [2].  At  the  other  extreme,  the  ion  and 
plasma  propulsion  devices  can  have  a  very  high  specific  impulse,  but  they  must  carry  a  heavy 
electrical  energy  source  to  deliver  the  power  necessary  for  high  ejection  velocities,  which  results 
in  a  low  thrust  value  [3].  Thus,  the  specific  impulse  of  the  advanced  electrical  rocket  engines 
can  be  appreciably  higher  than  those  of  chemical  or  some  nuclear  fission  rocket  engines.  This 
means  that  electrical  rocket  engines  need  to  carry  relatively  little  propellant  because  the 
propellant  mass  is  ejected  at  a  very  high  velocity.  The  low  thrust  values  of  the  electrical  rocket 
engines  imply  that  they  are  not  useful  in  the  fields  of  strong  gravitational  gradients  (such  as  for 
take-offs  or  landings  on  earth)  [5J. 

The  new  U.S.  space  mission  to  explore  Mars  initiating  from  the  Earth  orbit  would  require  a 
1400-sec  Isp  engine  and  a  longer  duration  of  flight  [6].  Naturally,  new  concepts  beyond  the 
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chemical  and  other  conventions!  propulsion  scenarios  discussed  above  would  be  required  to 
satisfy  the  higher  Lp  and  the  sufficient  flight  duration. 

Advanced  concepts  for  non-conventional  propulsion  developments  which  include  fusion 
energy,  anti-matter  energy,  etc.  have  been  proposed  [7].  One  of  the  most  promising  propulsion 
systems,  which  minimizes  the  size  of  the  device  and  the  total  mass,  is  considered  to  be  a  fusion 
rocket.  The  attractive  features  of  a  fusion  rocket  are  its  capabilities  to  convert  the  kinetic 
energies  of  energetic  charged  particles  directly  to  the  thrust  power  for  propulsion.  As  illustrated 
in  Table  2  and  Fig.  1,  the  yield  per  unit  mass  from  chemical  and  nuclear  fission  are  orders  of 
magnitude  smaller  than  those  of  nuclear  fusion.  Typical  specific  impulse  for  fusion  rockets 
could  range  up  to  10,000  sec  while  maintaining  relatively  high  thrust  levels. 


Table  2.  Yield  From  Various  Energy  Sources 

I 

ENERGY  RELEASE(JZkg) 

CHEMICAL 

Conventional  ( LO2/LH2 ) 
Recombination  (H  +  H* H2) 

Water,  Hydrogen 

Hydrogen 

1.5xl07 

2.18x10® 

NUCLEAR  FISSION 

U^.U^.Pu239 
(~200MeV  /  U235  fission) 

Radioactive 

Fission  Fragments 

Neutrons,  Gammas 

8.2xl013 

NUCLEAR  FUSION 

DT  ( 0.4/0.6 ) 

Helium,  Neutrons 

3.38xl014 

CatDD(l.O) 

Hydrogen,  Helium  & 
Neutrons 

3.45xl014 

DHe3  (0.4/0.6) 

Hydrogen,  Helium 
( Very  Few  Neutrons ) 

3.52xl014 

Recently,  the  Phillips  Laboratory  completed  a  fusion  propulsion  study  [8]  in  which  a  50-m 
long  translating  compact  torus  (TCT)  (e.g.,  moving  spheromak  reactor),  operating  with  D-3He 
fusion  fuel,  was  shown  to  be  a  prime  fusion  rocket  candidate.  The  large  size  of  the  TCT  device, 
delivering  a  total  payload  mass  of  36,000  kg  for  the  round  trip  between  a  low  Earth  orbit  (LEO) 
and  the  Geosynchronous  Earth  orbit  (GEO),  implied  that  the  thrust-to-weight  ratio  (i.e., 
condition  for  proper  acceleration)  would  be  rather  low.  The  ignition  condition  for  D-3He  fuel 
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PROPULSION  SYSTEM  PERFORMANCE  (LATER) 


Figure  1.  Propulsion  System  Specific  Mass  (kg/kWj)  [29] 


fusion  propulsion  device.  NASA  has  alrorecnulyrecoiranrndwl  that  a  D-3He  fueled  field- 
icversed  configuration  (FRO ,  another  compact  locus  device,  be  their  duke  for  fusion  space 
propulsion  since  the  FRC  could  provide  both  desired  characteristics  of  high  beta  (  90  %)  (Lc^ 
the  ratio  of  plasma  pressure  to  magnetic  field  pressure)  and  good  conversion  of  energy  to  thrust 
[9].  However,  the  DOE’s  FRC  experiment  is  no  longer  supported  and  further  development  of 
this  concept  would  be  somewhat  hampered  if  NASA  decides  to  pursue  it 

One  of  the  most  interesting  and  practical  plasma  thruster  concepts  is  the  magneto-plasma- 

dynamic  (MFD)  are  in  which  plasma  is  being  driven  by  magnetic  fields  generated  by  the 

currents  in  the  arc  [5].  The  thruster  geometry  is  coaxial,  with  a  cylindrical  cathode  at  the  center 

and  an  annular  anode  around  it  The  two  electrodes  are  separated  at  the  back  of  the  device  by  an 

insulating  plate  as  illustrated  in  Fig.  2.  The  current  from  the  ring  anode  to  die  cathode  produces 

an  azimuthal  magnetic  field  which  exerts  a  pressure  against  the  arc  plasma.  Propellant  from  gas 

— >  — > 

inlet  flows  through  die  arc  where  it  is  ionized  and  faced  away  by  die  magnetic  field  via  a  J  X  8 
interaction.  The  ampler  illustration  of  the  basic  principle  of  die  MPD  effect  is  illustrated  at  the 
lower  portion  of  Fig.2.  Efficiencies  as  high  as  50%  [10]  and  higher  specific  impulse  values  up  to 
10,000  sec  were  obtained  at  low  values  of  propellant  flow.  Current  MPD  tests  at  the  Phillips 
Laboratory  operate  in  a  2  msec-pulse  mode  with  die  capacitor  bank  charged  to  800  volts  and 
delivering  40,000  amps.  Instruments  are  used  to  measure  the  current  and  voltage  seen  by  the 
MPD  thruster  during  each  firing.  Typical  range  of  plasma  temperature  in  MPD  is  a  few  eV. 
However,  in  order  to  create  fusion  conditions  with  D-3He  fuels  in  an  MPD-type  coaxial  device, 
much  higher  current  and  plasma  temperature  (in  keV  range)  are  needed  which  necessitate  a 
dense  plasma  focus  device. 

The  dense  plasma  focus  (DPF),  one  of  the  complementary  fusion  devices  between  low-beta 
tokamak  and  high  compression  inertial  fusion,  is  a  system  of  coaxial  electrodes  which  allows  the 
formation  and  subsequent  propagation  of  a  thin  plasma  sheath  in  the  annular  region  between  the 
center  anode  and  the  outer  cathode.  When  the  sheath  reaches  the  end  of  the  anode,  magneto- 
hydrodynamic  (MHD)  instabilities  develop  and  the  sheath  disconnects  from  the  cathode  and 
collapses  toward  the  axis  forming  a  very  small  region  of  high  density,  hot  plasma  like  in  a  linear 
pinch  device  [11].  It  is  here  where  fusion  reactions  take  place  and  generate  the  energy  to  be 
used  for  propulsion. 

The  DPF  is  similar  in  geometry  to  the  magnetoplasmadynamic  (MPD)  thruster  [5]  currently 
used  in  electric  propulsion  but  differs  greatly  in  its  operation.  There  are  two  common  plasma 
focus  types:  The  Mather  type,  like  a  coaxial  plasma  gun,  has  the  current-carrying  plasma  sheath 
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Figure  2.  The  MPD  Thruster 


accelerated  axially  [12],  and  die  Filippov  type  with  the  plasma  sheath  accelerated  radially 
inward  [13].  The  DPF  typically  will  use  currents  which  are  about  1,000  times  greater,  and 
unlike  the  MPD  does  not  currently  operate  in  the  steady-state  mode.  The  MPD  forms  a  stable 
sheath  near  die  end  of  the  electrode  system  and  makes  no  use  of  a  rundown  phase  which  occurs 
in  the  DPF.  However,  the  key  distinction  between  the  two  devices  is  that  die  MPD  makes  no 
attempt  to  use  die  tremendous  amounts  of  energy  available  from  fusion.  It  is  here  that  the  DPF 
gains  an  advantage  over  most  pron»,oion  concepts.  For  example,  the  fusion  of  deuterium  (D) 
and  helium-3  ( 3  He)  can  release  almost  five  times  more  energy  per  unit  mass  than  the  fissioning 
of  a  uranium-235  nucleus. 

In  the  DPF,  the  plasma  sheath  is  initially  created  wher  a  large  current  is  discharged  through 
the  center  anode.  The  resulting  potential  difference  causes  the  current  to  arc  between  the 
electrodes.  In  the  process,  the  fill  gas  (fusion  fuel)  is  ionized  and  forms  an  azimuthally 
symmetric  plasma  sheath  in  the  annular  region  between  the  electrodes.  The  current  flowing 
through  the  ancde  also  produces  an  azimuthal  magnetic  field,  Be,  which  interacts  with  the 
plasma  sheath  current.  This  results  in  the  propagation  of  the  sheath  down  the  length  of  the  anode 
due  to  the  Jr  x  Be  force.  Figure  3  shows  the  cylindrical  thruster  configuration  as  well  as  the 
directions  of  the  current,  magnetic  field,  and  sheath  propagation.  During  "rundown,"  some 
fraction  of  the  fill  gases  is  entrained  in  the  sheath  and  carried  down  to  the  end  of  the  anode.  As 
the  sheath  reaches  the  end  of  the  anode,  it  does  collapse  or  "focus"  radially  inward  toward  the 
central  axis  of  the  device,  forming  a  high  density  (=  1026  m-3),  hot  plasma  where  fusion 
reactions  may  take  place.  This  number  density  may  change  depending  on  the  fraction  of  initial 
fuel  which  is  trapped  in  the  pinch  region. 

The  pinched  plasma  expands  and  contracts  several  times  before  it  eventually  becomes 
completely  disrupted  by  plasma  instabilities.  It  is  particularly  susceptible  to  the  m  =  0  "sausage" 
and  m  =  1  "kink"  instabilities.  Figure  4  gives  a  graphic  representation  of  the  different  phases 
which  occur  during  one  cycle  of  the  DPF.  The  pinch  lifetime  is  typically  very  short,  on  the  order 
of  a  microsecond  [14].  However,  if  the  pinch  lifetime  can  be  made  sufficiently  long  to  allow  a 
good  fusion  bum  inside  the  pinch,  the  DPF  could  provide  enough  energy  to  propel  spacecraft  at 
either  high  thrust,  high  specific  impulse,  or  both.  While  the  rundown  can  be  predicted  with 
reasonable  accuracy,  the  collapse  and  subsequent  plasma  behavior  are  not  well  understood  and 
are  in  great  need  of  further  study. 
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1)  BREAKDOWN 


2)  RUNDOWN 


3)  FOCUS  (PINCH) 


1)  RE-BREAKDOWN 


Figure  4.  Phase  of  the  DPF  Cycle 


There  have  been  numerous  activities  with  the  plasma  focus  experiments  around  the  world 
[15-17].  A  break  even  plasma  focus  reactor  concept  has  been  proposed  by  Herald  and  Hayd 
operating  on  3.0  MJ  energy  at  1.8  MV  voltage  with  6.0  MA  current  The  Stuttgart  group 
proposed  a  power  level  range  of  lOOkW  to  megawatts  for  fusion  reactions  and  plasma  stability 
studies  [15]. 

In  previous  experiments  at  the  Lawrence  Livermore  National  Laboratory,  various  mixtures 
of  D2  and 3  He  were  used  with  two  plasma-focus  devices  that  operated  with  a  stored  energy  of  73 
kJ  at  18  kV,  and  with  67  kJ  at  70  kV  [18].  Though  the  purpose  of  the  experiment  was  to  develop 
a  diagnostic  tool  for  determining  the  plasma  ion  temperature  by  evaluating  the  D-D/  D-3He 
yield  ratio,  the  significant  yield  of  14.6  MeV  protons  from  the  3He  (d,p)4He  reaction  was 
observed.  For  a  50/50  D2/3He  mixture,  the  D-D/  D-3He  yield  ratio  corresponded  to  deuterium 
beams  with  28  to  47  keV  incident  upon  stationary  3  He  ions. 

Also  a  1-MJ  plasma  focus  device  was  used  at  Frascati  to  measure  the  D-D/  D-3He  yield 
ratio  [19].  A  new  calibration  of  the  (p,n)  cross  section  and  the  ratio  of  fusion  reaction  yields 
resulted  in  the  higher  deuteron  beam  temperature  at  72  keV.  Filippov  has  also  observed  high 
energy  ions  emitted  from  his  plasma  focus  device  [13];  Gullickson  and  Sahlin  have  measured 
more  than  1015  deuterons  above  330  keV  and  1012  deuterons  with  energy  greater  than  5  MeV 
[20];  and  Mozer  et  al.  have  reported  the  fast  deuterons  with  energies  greater  than  350  keV  by 
using  50  kJ  at  18  kV  plasma  focus  device  [21].  Quasisteady  multimegawatt  MPD  thruster 
performance  and  some  comparative  analaysis  of  large  plasma  focus  experiments  (of  360  kJ  and 
500  kJ  devices)  performed  at  IPF  (Stuttgart)  and  at  EPJ  (Poland)  were  also  reported  [22-25]. 

More  specific  experimental  component  studies  with  a  chamber  magnetic  nozzle  provided  30% 
increase  of  the  ion  density  with  8%  addition  of  the  electric  power  to  the  discharge  chamber  [26]. 
Small  discharge  chamber  length  yielded  high  extracted  ion  fractions  [27]  and  the 
characterization  of  plasma  flow  through  magnetic  nozzles  was  also  performed  by  Gerwin,  et  al. 
[28]. 

Various  plasma  focus  experimental  data  were  compiled  in  an  effort  to  understand  the  scaling 
relations  between  neutron  yield  vs.  pinch  current  and  stored  energy,  and  they  are  illustrated  in 
Fig.  5  [30].  It  is  apparent  from  the  above  observations  that  the  plasma  focus  device  is  capable  of 
igniting  D-D  and  D-3He  fuels.  Hence,  the  goal  of  this  study  was  not  only  to  recognize  the 
MPD-type  plasma  focus  device  but  also  to  optimize  engineering  parameters  including  optimum 
capacitor  bank  voltage,  energy,  inductance,  and  plasma  temperatures.  The  optimized  parameters 
will  allow  the  maximum  energy  transfer  from  the  capacitor  bank  to  the  plasma  pinch. 
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Figure  5.  Neutron  Yield  as  Functions  of  Pinch  Current  and  Stored  Energy  [Ref.  30] 


Plasma  focus  calculations  in  the  literature  were  also  compared.  Typical  calculational  models 
include  1)  a  simple  LRC  circuit  model  in  which  entire  plasma  focus  is  included  in  the  circuit 
equation  by  regarding  it  as  contributing  a  variable  inductance;  2)  a  plasma  focus  model  in  which 
the  plasma  focus  is  divided  into  three  areas  of  i)  radial  lift-off  region  above  insulator,  where 
current  is  generated  due  to  a  pressure  balance  between  magnetic  field  and  particle  pressure,  ii) 
run-down  region  above  anode  (with  snow-plow  model),  and  iii)  pinch  region  off  the  end  of  the 
anode;  and  3)  a  2-D  MHD  model  in  which  both  energy  and  momentum  are  conserved.  The  2-D 
MHD  model  was  developed  [31]  and  used  by  Eltgroth  [32]  to  check  experimental  values  of  the 
Livermore-I  &  -II  and  Frascati  devices.  Other  related  numerical  studies  were  done  for  a  soft  x- 
ray  experiment,  high  temperature  pinch,  current  distribution,  and  MPD  thrusters  [33-38]. 

Various  scaling  relations  were  also  studied  by  others  using  the  Boltzmann  equation,  Maxwell 
equations  and  binary  collision  approximation  [39].  However,  the  aspects  of  plasma  scaling 
relations  will  not  be  included  in  this  current  report. 

In  general,  all  the  calculational  models  are  interconnected  with  the  conservation  of  energy 
and  momentum.  2-D  MHD  models  are  more  rigorous  and  often  rely  on  complex  numerical 
codes.  However,  by  developing  a  simple  plasma  focus  model  with,  for  example,  a  plasma  snow 
plow  driven  by  magnetic  pressure  during  the  run-down  phase  [32],  one  can  simplify  the 
computations  analytically.  The  present  study  employees  the  analytical  plasma  focus  model  using 
the  steady-state  MHD  momentum  equation  coupled  with  the  equations  for  maximum  attainable 
plasma  current  described  in  terms  of  capacitance,  inductance,  charging  voltage,  etc. 
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H.  FEASIBILITY  STUDIES  OF  DENSE  PLASMA  FOCUS  AS  A  PROPULSION 

THRUSTER 

A  systems  model  of  DPF  is  developed  for  both  conventional  and  spin-polarized  D-3He  fuels 
to  investigate  the  feasibilities  of  DPF  as  a  propulsion  thruster.  Three  modes  of  operations  (i.e., 
two  pulsed  operations  with  and  without  hydrogen  propellant,  and  one  impulsive  firing  mode  with 
hydrogen  propellant)  were  studied  and  each  was  investigated  for  its  usefulness  in  space  travel 
with  special  attention  paid  to  a  manned  Mars  mission.  This  portion  of  the  study  with  the 
conventional  D-3He  is  described  in  Section  II.  1.  The  spin-polarized  D-3He  operation  was  also 
analyzed  to  study  the  effect  of  polarized  fuels  on  the  systems  operations  on,  for  example,  thrust- 
to-weight  ratio  (F/W)  and  specific  impulse  (Isp).  Increased  fusion  power  and  decreased  radiation 
losses  for  the  spin-polarized  case  provided  the  increased  values  of  F/W  and  Isp  which  are 
described  in  Section  II. 2. 

H.  1.  PARAMETRIC  STUDIES  OF  DENSE  PLASMA  FOCUS  WITH  D-3He  FUEL* 

H.1.1  DPF  COMBUSTOR  MODEL 

The  rundown  velocity,  U^,  can  be  predicted  accurately  by  solving  the  steady-state 
••  momentum  equation  for  the  plasma  sheath  neglecting  dissipative  effects,  starting  with  [32] 

pi^(0„0„)— vp  «) 


then  taking  only  z-components  and  integrating  gives 


U, 


run 


I  Ho 
\  8*2  r»Pi 


(2) 


where  I  is  the  current  discharged  through  the  anode,  ra  is  the  radius  of  the  anode  and  Pi  is  the 


*  An  earlier  version  of  this  section  of  the  report  has  been  published  as  PL-TR-91-3014  and  the  revised 
version  has  been  documented  as  the  M.S.  project  report  by  C.  Leakeas  as  a  partial  requirement  for  the 
M.S.  degree  from  Purdue  University. 
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initial  fill  gas  density.  The  sheath  then  teaches  the  end  of  the  anode  and  collapses  forming  a 
small,  hot  plasma.  (It  will  be  shown  later  that  the  final  temperature  depends  on  many  factors 
including  the  capacitor  discharge  current.  Also,  the  pinch  dimensions  are  assumed  to  be 
independent  of  operating  conditions.)  If  one  assumes  that  a  fraction,  f,  of  particles  goes  into  the 
pinch,  and  makes  a  rough  estimate  of  the  dimensions  of  the  resulting  pinch  formation,  one  can 
determine  the  number  density  of  particles  inside  the  pinch.  To  determine  the  temperature  inside 
the  pinch,  one  assumes  a  balance  between  plasma  pressure  and  the  magnetic  pressure  due  to  the 
external  azimuthal  field,  Be,  generated  by  the  current  in  the  pinch,  as  in  Eq.  3. 


npkT  = 


-5L 

2Mo  ’ 


(3) 


where  kT  is  the  product  of  Boltzmann's  constant  and  the  plasma  temperature  (in  degrees)  and  p<, 
is  the  permeability  of  free  space. 

By  Ampere’s  Law,  Be  at  the  pinch  surface  is 


Be  * 


Mol 

2nrp 


(4) 


Solving  for  the  pinch  number  density,  np,  in  terms  of  the  initial  fill  density  and  pinch  and 
electrode  dimensions,  one  finds 


fpjl,(rf-r*) 
lp  rp  ntp 


(5) 


where  f  is  the  fraction  of  initial  fuel  which  is  trapped  in  the  pinch  region,  1,  and  lp  are  the  anode 
and  pinch  lengths,  ra,  rc,  and  rp  are  the  anode,  cathode,  and  pinch  radii,  respectively,  and  nip  is 
the  average  mass  of  particles  ;<i  the  pinch.  The  fraction,  f,  is  left  as  an  independent  variable 
dunng  thruster  evaluation,  but  is  assumed  to  be  about  17.5%  for  the  baseline  case.  This  value  is 
atr.ed  at  by  assuming  that  about  70%  of  all  fuel  is  entrained  during  rundown,  and  of  that,  25% 
is  captv  ed  in  the  pinch.  8  jbstituting  Eqns.  4  and  5  into  Eqn.  3  gives  an  expression  for  the 
plasma  temperature  inside  the  pinch. 
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i.nn  _  Mo  I2  Mp  lp 

““  ftfthl.fr* -r  \)' 


(6) 


This  gives  the  plasma  temperature  for  any  current  I.  The  maximum  attainable  current  must  now 
be  calculated  as  a  function  of  the  electrical  parameters  of  the  system  (e.g.  capacitance, 
inductance,  charging  voltage,  etc.)-  The  maximum  attainable  current  is  given  by 


where 


In™  =0.64 


(7) 


(8) 


(9) 


K  J 


(10) 


where  C  is  the  total  capacitance,  V  is  the  charging  potential,  and  1*  is  the  initial  circuit 
inductance  [40].  Thus,  from  the  initial  conditions  defined  in  Eqs.  8-10,  the  maximum  current 
and  resulting  plasma  pinch  temperature  can  be  found  with  Eqs.  6  and  7.  It  is  assumed  through 
the  rest  of  this  report  that  the  plasma  pinch  temperature  continues  to  scale  as  current  squaied  as 
shown  in  Eqn.  6  and  Figure  9,  although  this  scaling  seems  to  fail  for  currents  above  1 MA  due  to 
saturation  and  degradation  effects  [25]. 

The  plasma  focus  device  analyzed  in  this  report  is  assumed  to  be  identical  to  the 
"Livermore-I"  dense  plasma  focus.  Therefore  the  same  geometrical  and  electrical  parameters 
which  were  used  in  the  operation  of  this  device  will  be  adopted  [20,32].  The  parameters  used 
can  be  found  in  Table  3  where  asterisks  denote  assumed  values.  If  operated  at  the  values  in 
Table  3,  the  Livermore-I  focus  should  be  capable  of  a  maximum  current  of  1.245  MA  and  a 
maximum  plasma  temperature  of  about  300  eV.  This  plasma  temperature  is  much  too  low  to 
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produce  significant  amounts  of  fusion  energy.  As  seen  in  Figure  6,  one  would  ideally  operate  at 
kT  greater  than  about  SO  keV  (depending  on  the  fuel)  in  order  to  maximize  the  reaction  rate 
parameter.  Therefore  the  capacitor  banks  are  assumed  to  be  capable  of  delivering  in  excess  of 
20  MA.  Using  the  values  in  Table  3,  the  DPF  performance  can  be  modeled  for  a  wide  range, of 
currents  using  the  simple  scaling  laws  found  in  Eqs.  4-6. 


Table  3.  Livermore*!  Dense  Plasma  Focus  Parameters 

V-  27,000  V 

C  *  3.55  x  10"4  F 

Lo  =  2.5  x  10-8  H 

1,  *  0.382  m 

r,  *  0.0508  m 

rc  =  0.080  m 

*  f«0.175 

*lp=  0.0254  m 

*  rp  ■  0.0015  m 

pi «  2.2  x  10T4  kg/m3 

♦Denotes  assumed  values. 

In  order  to  complete  the  development  of  the  DPF  model,  the  fusion  fuels  need  to  be  defined. 
There  are  several  desirable  criteria  to  be  considered  in  the  choice  of  fusion  fuels.  The  fusion 
fuels  should  have  a  high  power  density.  Of  all  the  fusion  reactions  being  considered  today,  the 
D-3He  reaction  has  the  highest  power  density  at  18.3  MeV  per  reaction.  D-T  is  second  at  17.6 
MeV, 3  He-3  He  releases  12.9  MeV,  and  p-u  B  releases  8.7  MeV  per  reaction. 

The  fuel  should  have  a  high  reaction  rate  parameter  at  low  temperature  to  achieve  good 
fusion  bum  before  radiative  losses  overwhelm  the  system.  The  reaction  rate  parameter,  <ov>,  is 
a  function  of  plasma  temperature  and  determines  how  quickly  these  reactions  proceed  at  a  given 
temperature.  Figure  6  shows  the  reaction  rate  parameters  for  some  typical  fusion  fuels  [41],  D- 
T  has  the  highest  reaction  rate  parameter  at  low  temperatures  and  the  D-3He  reaction  rate 
parameter  is  slightly  less.  The  other  "advanced"  fuels  such  as  p-6Li,  p-n  B,  and  3 He-3 He  must 
be  operated  at  very  high  temperatures  [42,43]  and  may  be  impractical  because  synchrotron 
radiation  losses  increase  in  proportion  to  the  electron  temperature  squared  at  the  lowest  order, 
thus  making  ignition  difficult.  The  D-D  reaction  has  a  reaction  rate  parameter  which  is  slightly 
less  than  that  for  D-T,  but  this  fuel  has  two  major  disadvantages.  It  has  a  much  lower  power 
density  (about  4  MeV)  and  has  50%  probability  of  producing  a  neutron  with  each  reaction. 
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Figure  6.  Reaction  Rate  Parameters  for  Various  Fuels 
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Since  neutrons  cannot  be  directed  to  produce  thrust  by  a  magnetic  field,  it  is  desirable  to 
minimize  their  production.  An  ideal  reaction  would  release  all  of  its  energy  in  the  form  charged 
particles.  Of  the  "easily"  ignitable  fuels,  D-T  releases  about  80%  of  its  energy  as  neutrons,  while 
the  D-3He  reaction  releases  no  neutrons.  However,  the  secondary  background  D-D  in  D-3He 
reaction,  which  releases  about  75%  of  its  energy  in  neutron  radiation,  can  contribute  significant 
neutron  production  to  the  D-3He  reaction.  The  other  "advanced"  fuels  mentioned  earlier  are 
considered  to  be  completely  a  neutronic. 

After  considering  each  fuel’s  characteristics  in  these  three  areas,  D-3He  was  chosen  as  the 
fuel  to  be  used  in  this  DPF  study.  Its  leaner  neutron  production  (compared  to  D-T),  high  power 
density  and  high  reaction  rate  parameter  at  relatively  low  temperatures  were  key  factors  in  its 
selection. 

n.1.2  DPF  PROPULSION  SYSTEM  MODEL 

The  DPF  fusion  propulsion  system  (see  Figure  7)  consists  of  the  feed  and  cooling  system,  the 
electrical  power  system  and  the  thruster  system.  The  feed  and  cooling  system  consists  of  three 
tanks  for  the  hydrogen,  helium-3,  and  deuterium,  associated  plumbing  to  control  and  direct  the 
flow  of  these  gasses,  and  the  associated  coolant  passages.  Also  necessary  in  this  system  are  a 
number  of  pumps  to  drive  propellant  and  fuel  fiow,  as  well  as  to  provide  the  pressure  necessary 
for  the  coolant  to  enter  the  high  pressure  side  of  the  turbine.  The  deuterium  and  helium-3  are 
used  as  the  fuel  to  drive  the  thruster.  The  hydrogen  is  used  for  cooling,  driving  the  turbine,  and 
may  then  also  be  used  as  propellant  to  provide  increased  thrust. 

The  electrical  system  consists  of  a  turbine,  electric  generator  and  the  capacitor  banks 
necessary  to  produce  the  large  current  pulses  which  are  required  by  the  thruster.  The  electricity 
produced  by  the  turbo-generator  is  used  to  meet  system  requirements  and  to  help  recharge  the 
capacitor  banks  for  each  shot. 

The  thruster  system  consists  of  the  DPF  itself  as  well  as  a  mixing  chamber  and  a  magnetic 
nozzle  if  one  chooses  to  operate  the  DPF  at  very  high  propellant  temperatures.  As  will  be  seen 
later,  the  magnets  necessary  for  our  purpose  are  relatively  small  and  constitute  a  small  fraction 
of  total  system  mass.  The  mixing  chamber  is  only  necessary  if  the  DPF  propulsion  system  is  to 
be  operated  with  hydrogen  propellant.  It  is  a  hollow  cylindrical  cavity  where  the  fusion  reaction 
products  will  mix  with  cold  hydrogen  propellant.  It  is  assumed  that  the  resulting  mixture  leaves 
the  chamber  with  a  uniform  temperature  and  produces  thrust. 
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PLASMA  FOCUS  FUSION 
PROPULSION  SYSTEM 


Figure 


7.  Plasma  Focus  Fusion  Propulsion  System 


The  exhaust  power  produced  in  any  type  of  rocket  engine  is  given  by 


Fra  —  T  n^propeDjni 


(U) 


v/heic  F  is  the  thrust  and  Ura  is  the  propellant  exhaust  velocity.  If  one  assumes  that  the  thrust  is 
parallel  to  die  exhaust  velocity  and  defines  die  specific  impulse  by 


I*  = 


U 


ex 


g 


(12) 


and  the  thrust  by 


P  =  Olpropdknt  De*  » 


(13) 


the  exhaust  power  can  then  be  written  as 


— 2  g  ^  Vp 


(14) 


where  g  is  the  gravitational  acceleration  and  specific  impulse  is  a  measure  of  how  efficiently 
propellant  is  used  [44].  Eqn.  14  shows  the  competing  nature  of  thrust  and  specific  impulse. 

Both  are  desirable,  but  for  a  fixed  engine  power  an  increase  in  one  requires  a  decrease  in  the 
other.  However,  with  the  high  exhaust  powers  attainable  with  fusion,  it  should  be  possible  to 
attain  reasonably  high  values  of  both  parameters  simultaneously. 

Three  possible  modes  of  operation  for  the  DPF  propulsion  system  were  investigated: 

1 )  Pulsed  operation  of  the  DPF  for  long  periods  of  time  with  no  hydrogen  propellant 
exhausted.  The  fusion  products  are  produced  and  immediately  expelled  to  produce  thrust. 
The  total  time  that  the  thruster  is  fired  is  comparable  to  the  total  trip  time. 

2)  Pulsed  operation  of  the  DPF  for  long  periods  of  time  with  the  addition  of  moderate 
quantities  of  hydrogen  propellant.  The  hydrogen  is  used  to  provide  electric  power  and 
also  provides  additional  thrust  because  of  increased  mass  flow  rate  in  the  exhaust  with 
some  loss  of  Isp. 
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3)  Pulsed  operation  of  the  DPF  for  short  periods  of  time  during  which  large  quantities  of 
hydrogen  are  exhausted  in  a  high  thrust  impulsive  bum.  This  "impulsive"  bum  reduces 
gravitational  losses,  makes  much  higher  thrust-to- weight  ratios  possible  and  would  most 
likely  be  used  during  interplanetary  travel. 

Mode  1)  Pulsed  operation  with  no  hydrogen  propellant 

This  mode  involves  a  closed  coolant  cycle  and  would  therefore  require  large  radiators  (0.07 
kg/kWe)  [45]  to  dissipate  the  heat  produced  by  resistive  heating  of  the  electrodes  and  radiative 
power  losses  in  the  thruster  walls.  It  is  still  possible  to  generate  electricity  from  the  turbo¬ 
generator  before  the  coolant  enters  the  radiators.  The  calculations  made  for  this  mode  are 
identical  to  those  made  using  the  FORTRAN  code  in  appendix  of  this  report,  except  that  one 
must  include  the  radiator  mass  and  use  a  zero  propellant  mass  flow  rate.  The  system  is  similar  to 
that  shown  in  Figure  7,  but  with  the  addition  of  radiators  in  the  coolant  loop  and  removal  of  the 
mixing  chamber. 

When  the  fusion  fuels  react  in  the  pinch,  it  is  assumed  that  very  little  charged  particle  power 

is  retained  in  the  pinch.  Thus,  when  the  charged  particles  leave  the  pinch,  their  energies  are 

v  * 

known  simply  as  a  function  of  the  fusion  fuels  used  (D-  He  produces  14.7  MeV  protons  and  3.6 
MeV  alpha  particles).  The  velocities  of  these  particles  are  very  high  (some  over  107  m/s). 

These  velocities  lead  to  specific  impulse  values  on  the  order  of  106  s.  However,  because  of  the 
low  mass  flow  rate  exiting  the  pinch,  thrust  values  are  on  the  order  of  about  44.5  N  (10  lbf),  and 
the  main  contribution  to  this  thrust  comes  from  the  expulsion  of  fill  gases  which  are  not  trapped 
during  the  rundown  and  pinch  phase  of  operation.  For  a  manned  Mars  mission  with  a  payload 
dominating  mass  of  105  kg  [46],  the  system  thrust-to- weight  ratio  (F/W)  upper  bound  is  about 
5.0  x  10"5.  If  the  addditional  mass  of  radiators,  shielding,  capacitors,  tanks,  fuel,  etc.  are 
considered,  F/W  decreases  further.  This  F/W  value  is  many  orders  of  magnitude  less  than 
conventional  chemical  rockets.  In  this  mode,  the  DPF  is  comparable  in  performance  to  electric 
propulsion.  Although  these  thrust  levels  have  applications  to  certain  types  of  missions  (perhaps 
orbital  transfer),  manned  interplanetary  travel  requires  larger  mission  Av’s  and  shorter  trip  times 
to  reduce  exposure  to  cosmic  radiation  and  weightlessness.  Therefore  this  mode  was  not 
considered  beyond  the  conceptual  state. 

Mode  2)  Pulsed  operation  with  hydrogen  propellant 

One  way  to  increase  F/W  values  is  by  exhausting  the  heated  coolant  to  increase  the  mass 
flow  rate  and  corresponding  thrust  given  in  Eq.  1 3.  In  doing  this  one  accepts  the  penalty  of 
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decreased  specific  impulse  as  a  necessary  means  to  increase  thrust.  The  increase  in  thrust  will 
allow  the  DPF’s  use  in  a  wide  variety  of  missions,  whereas  operating  without  hydrogen 
propellant  restricts  the  type  of  missions  for  which  it  can  be  used. 

The  Mode  2  system  is  illustrated  schematically  in  Figure  7.  As  the  capacitor  banks  are  being 
discharged  through  the  center  anode,  the  fuel  is  injected  and  is  caught  up  in  the  plasma  sheath’s 
rundown  as  illustrated  in  Figure  4.  The  plasma  collapses  and  pinches  at  the  end  of  he  anode  and 
produces  large  amounts  of  charged  fusion  products  as  well  as  Bremsstrahlung  and  synchrotron 
radiation  and  neutrons.  It  was  assumed  in  this  study  that  the  Bremsstrahlung  radiation,  which  is 
emitted  in  the  UV  spectrum,  is  completely  lost  However,  part  of  the  synchrotron  radiation  was 
assumed  to  be  absorbed  by  the  plasma  and  part  absorbed  in  the  walls  of  the  electrodes  and 
mixing  chamber.  The  portion  absorbed  in  the  walls  and  electrodes  was  assumed  to  be  about 
20%  for  the  baseline  case,  but  was  left  as  a  variable  parameter  to  observe  its  effects  on  thruster 
performance.  The  heat  generated  by  synchrotron  radiation  and  ohmic  heating  is  then  cooled  by 
the  flow  of  cold  liquid  hydrogen  propellant.  Because  of  material  limitations,  the  turbine 
entrance  temperature  was  constrained  to  be  no  greater  than  2,000  °K.  This  would  require 
advances  in  material  sciences  since  current  materials  restrict  temperatures  to  less  than  about 
1,100  °K  [47].  This  inlet  temperature  constraint  then  fixes  the  minimum  mass  flow  rate  of 
coolant  which  enters  the  mixing  chamber.  The  gas  is  then  expanded  through  a  turbine  used  to 
run  a  generator  which  recharges  the  capacitor  bank.  Complete  recharging  of  the  capacitors  is 
only  possible  at  higher  powers  and  larger  coolant  flow  rates.  The  flow  from  the  turbine  is  then 
used  as  propellant  to  absorb  the  energy  of  the  charged  particles  produced  in  the  pinch.  An  open 
cycle  was  chosen  to  avoid  heavy  radiators  that  would  greatly  increase  total  system  mass.  The 
propellant  was  assumed  to  absorb  all  of  the  fusion  product  energy  after  accounting  for  radiative 
losses  and  become  completely  dissociated  and  ionized  in  the  mixing  chamber.  The  resulting 
propellant  plasma  was  assumed  to  come  to  a  uniform  temperature  before  it  enters  a  meridional 
magnetic  nozzle  (axial  field  only).  The  magnetic  nozzle  (see  Figure  8),  which  would  require  a 
maximum  field  of  about  2  T  [28],  then  further  accelerates  the  particles  out  the  exit  to  even  higher 
velocities. 

In  doing  the  analysis,  several  simplifying  assumptions  were  made.  These  assumptions 
concerned  aspects  of  the  pinch  as  well  as  advancements  in  other  technologies  applicable  to  the 
DPF.  These  assumptions  are: 
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1)  Plasma  pinch  temperature  scales  as  current  squared. 

2)  Since  no  accurate  measurements  of  actual  pinch  dimensions  have  been  made,  a  rough 
estimate  was  used. 

3)  Ions  come  to  thermal  equilibrium  inside  the  pinch  allowing  the  use  of  Maxwellian  reaction 
rate  parameters. 

4)  Materials  will  be  developed  that  can  withstand  temperatures  much  higher  than  currently 
possible.  This  would  be  necessary  in  the  turbine  and  in  the  walls  of  the  mixing  chamber  to 
minimize  damages  due  to  high  heat  fluxes. 

5)  Electrodes  and  mixing  chamber  walls  can  be  sufficiently  cooled  to  prevent  damage.  Film 
cooling  may  be  possible,  but  at  the  cost  of  I,p. 

6)  Propellant  becomes  completely  dissociated  and  ionized  in  the  mixing  chamber  at 
5,000  °K. 

7)  Advances  in  capacitor  bank  technology  will  increase  specific  energies  by  a  factor  of  10 
and  allow  for  discharging  rates  of  100  Hz.  Capacitors  based  on  present  technology  offer  a 
specific  masses  of  about  0.2  kJAcg  [8]. 

8)  Confinement  times  can  be  increased  about  a  hundred  times  (to  about  10*4  s)  to  allow  for  a 
good  fusion  bum  (around  40%).  Since  reaction  rates  are  determined  by  plasma 
temperature,  longer  confinement  times  allow  for  more  fuel  to  be  burned. 

9)  Any  magnetic  fields  applied  downstream  do  not  adversely  affect  the  pinch  formation  or 
confinement  time. 

With  these  assumptions,  thruster  performance  was  investigated  while  varying  current,  fraction  of 

particles  trapped  in  pinch,  capacitor  bank  specific  energy,  total  firing  time,  and  fraction  of 

synchrotron  radiation  absorbed  in  walls  and  electrodes. 


Baseline  case: 

f  *  Fraction  of  particles  trapped  in  pinch  0.175 

FRACT  =  %  synchrotron  radiation  absorbed  in  walls 

and  electrodes  0.20 

SPECEN  =  Capacitor  bank  specific  energy  2  0kJ/k? 

DAYS  =  Total  thruster  firing  time  30  days 
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For  the  cases  of  continuous  pulsed  operation  (Modes  1  and  2),  DAYS  is  defined  as  the  length 
of  time  which  the  thruster  is  fired.  The  baseline  case  assumes  that  the  thruster  is  fired  for  30 
days  and  that  this  is  comparable  to  total  trip  time. 

Figure  9  shows  the  pinch  temperature  dependence  on  current  for  several  values  of  f,  the 
fraction  of  particles  trapped  in  the  pinch.  Using  the  assumed  DPF  dimensions  and  initial  fill  g&s 
density,  f  «=  0.175  gave  pinch  number  densities  which  are  close  to  experimentally  determined 
values  (n  £1019  cm”3)  [14].  Although  lower  f  gives  lower  pinch  number  density,  Figure  9  shows 
lower  f  also  gives  higher  plasma  temperature.  Operating  at  the  very  high  temperatures  necessaary 
to  ignite  some  advanced  fuels,  such  as  3He-3He  and  p-nB,  may  not  be  feasible  since 
synchrotron  radiation  increases  as  T2. 

Figure  10  shows  the  resulting  propellant  temperature  at  the  entrance  to  the  magnetic  nozzles- 
as  a  function  of  current  for  various  values  of  f.  As  current  is  increased  past  a  certain  point,  the 
extra  fusion  power  produced  cannot  continue  to  raise  the  temperature  of  the  increased  coolant 
(and  therefore  propellant)  flow  which  must  be  supplied  due  to  greater  heat  flux  to  the  walls  and 
electrodes.  Therefore,  this  function  does  not  increase  indefinitely,  but  has  a  definite  maximum  at 
about  15  MA  for  baseline  values  of  f  and  FRACT.  This  corresponds  to  the  maximum  in  Specific 
impulse  for  the  baseline  case  in  Figure  11.  As  f  increases,  the  number  density  in  the  pinch 
increases,  so  the  resulting  pinch  temperature  decreases:  the  reaction  rate  decreases  because  of 
the  lowered  temperature,  and  there  is  less  fusion  output  to  beat  the  propellant.  A  balance  is 
established  between  plasma  and  magnetic  pressures,  so  for  higher  values  of  f,  and  therefore 
particle  density,  a  higher  current  is  required  to  bring  the  plasma  temperature  up  to  its  maximunri. 
As  the  maximum  current,  Imu,  increases,  the  propellant  mass  flow  must  correspondingly 
increase  to  cool  electrodes  of  the  focus  device.  Figure  12  shows  that  the  current  which 
maximizes  specific  impulse  (about  15  MA  for  the  baseline  case  in  Figure  11),  also  produces  a 
maximum  F/W  for  the  baseline  case.  As  current  is  increased  beyond  this  optimum,  capacitor 
mass  and  required  coolant  mass  increase  resulting  in  a  decrease  in  thrust-to- weight  ratio. 

System  F/W  ratios  are  calculated  by  the  program  in  the  appendix  taking  into  account  all  Systran 
masses.  As  seen  in  Figure  12,  vehicle  F/W  peaks  at  about  15  MA  and  reaches  almost  0.003  for 
the  baseline  case,  while  a  typical  value  for  a  manned  Mars  mission  using  an  impulsive  bum  is 
about  0.2  [46].  This  would  seem  to  be  the  optimum  operating  regime  for  the  DPF  in  this  mode 
operating  at  baseline  conditions,  as  it  maximizes  both  specific  impulse  and  F/W. 

Another  problem  is  in  the  area  of  capacitor  bank  technology.  Modem  capacitors  allow  a 
specific  energy  of  about  0.2  kJ/kg  [8].  However,  to  supply  the  necessary  currents  to  the  thruster 
and  magnet,  these  specific  energies  would  require  capacitor  masses  on  the  order  of  40,000  kg 
(about  40%  of  the  assumed  payload  fear  a  manned  Mars  mission).  Advancements  in  capacitor 
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Figure  9.  Plasma  Pinch  Temperature  vs.  Current 
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Figure  10.  Propellant  Temperature  vs.  Current 
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technology  might  allow  specific  energies  of  2.0  kJ/kg  and  would  make  the  total  systemmass 
much  smaller  and  allow  higher  thrust-to- weight  ratios  (see  Figure  13).  Further  increases  over 
2.0  kJ/kg  change  thrust-to- weight  ratios  only  slightly  because  at  high  specific  energies  system 
mass  is  dominated  by  propellant  and  payload. 

The  final  parameter  was  the  total  firing  time.  Figure  14  shows  the  expected  decrease  in  F/W 
as  firing  time  increases  because  of  the  increase  in  propellant  mass  which  must  be  carried.  Once 
again,  the  maximum  value  seems  to  occur  at  about  15  MA.  Although  thrust  is  increased,  the 
problem  with  continuous  operation  is  apparent.  In  this  high  1^,  low  F/W  mode,  mission  times 
become  extremely  long  and  thrust-to- weight  values  drop  even  further,  giving  the  DPF  limited 
usefulness  for  interplanetary  travel. 


Mode  3)  Impulsive  firing  with  hydrogen  propellant 

The  plasma  focus  propulsion  system  can  also  be  operated  by  firing  for  a  short  period  of  time 
while  exhausting  great  quantities  of  propellant  In  this  way,  the  propellant  has  been  exhausted 
and  is  no  longer  considered  to  contribute  to  the  total  system  mass.  This  decreased  system  mass 
allows  larger  acceleration  for  the  same  thrust  resulting  in  larger  F/W  ratios.  These  higher  F/W 
ratios  decrease  the  required  Av  for  a  given  mission  resulting  in  a  decrease  in  trip  time.  By 
adding  additional  hydrogen  flow  to  the  coolant  flow,  thrust  can  be  increased  considerably  but 
this  reduces  the  propellant  outlet  temperature  and  therefore  I^. 


A  key  issue  in  this  analysis  is  how  the  system  masses  are  obtained.  It  is  important  to  keep 
system  masses  low  to  improve  F/W  values  and  to  minimire  the  cost  of  raising  the  vehicle  from 
earth  into  low  earth  orbit  (LEO).  Payload  mass  is  a  constant  and  is  fixed  for  a  given  mission. 
Propellant  man  is  fixed  by  the  mission  Av  and  the  exhaust  velocity,  and  may  be  found  from  the 
rocket  equation: 


Mi  TF 
Mi 


(15) 


where  Mi  is  the  initial  mass  (i.e.,  total  system  mass,  payload  and  propellant)  and  Mf  is  the  final 
mass  after  the  bum  (total  system  mass  and  payload  only).  The  rocket  equation  can  be  simplified 
to 


*  g  I*p  In 


(16) 
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INITIAL  F/W  RATIO 


Figure  13.  Thrust-to-Weight  Ratio  vs.  Current 
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Because  of  the  possibility  for  high  the  plasma  focus  is  deariy  capable  of  high  Av’s  and  quick 

trip  tunes  if  adequate  thrast-to-weight  values  can  be  attained.  Figure  15  shows  the  total  Av 
capability  of  the  DPF  propulsion  system  as  a  function  of  payload  mass  fraction  at  1^  = 3,500  sec 
and  11 1,250  N  (25,000  lbf)  of  thrust  Another  very  importer  ystem  mass  is  the  capacitor 
mass.  It  is  necessaiy  to  project  die  capacitor  mass  by  scaling  from  some  reference  case.  The 
"LivermoreF  device  was  chosen  as  die  reference  to  remain  consistent  It  was  initially  assumed 
that  the  capacitor  mass  scaled  as  the  current  squared  The  results  presented  thus  far  and  to  be 
presented  here  have  been  obtained  with  the  I2  scaling.  However,  by  striving  Eqs.  7-10,  it  can  be 
shown  that  the  capacitor  m?ss  scales  as  Iw.  This  seemingly  small  difference  is  in  fact  very 
important  and  the  following  analysis  is  repeated  in  Appendix  A  using  Iw  scaling. 

In  this  impulsive  firing  mode,  it  is  important  to  find  the  optimum  propellant  mass  flow  rate  in 
order  to  maximize  F/W  while  maintaining  acceptable  specific  impulse.  Figures  16-20  show  the 
dependence  of  F/W  on  propellant  mass  flow  rate  for  various  currents  and  mission  velocity 
increments.  Figure  16  would  seem  to  show  that  F/W  can  be  increased  indefinitely  by  simply 
exhausting  more  propellant.  However,  this  lowers  the  propellant  outlet  temperature  which 
increases  the  total  amount  of  propellant  required  in  accordance  with  Eq.  IS.  This  is  not  apparent 
in  missions  with  small  Av’s  since  the  originally  required  propellant  mass  is  so  small  compared  to 
other  system  masses.  Figures  16  and  17  actually  show  cases  where  propellant  becomes  the 
dominant  mass  and  an  increase  in  propellant  mass  flow  rate  actually  decreases  system  F/W.  As 
a  result  of  the  l2  scaling  it  is  always  advantageous  to  run  the  thruster  at  the  highest  current 
possible  until  the  pinch  temperature  reaches  the  temperature  corresponding  to  maximum  in  the 
D-3He  reaction  rate  parameter  curve.  Optimum  F/W  ratios  for  the  high  Av  missions  are  still  an 
order  of  magnitude  below  the  desired  value  of  0.2  for  a  manned  Mars  mission.  Figure  21  shows 
the  variation  of  specific  impulse  for  Av  =  5  km/s,  but  is  in  fact  valid  for  any  mission  since  the 
propellant  exit  velocity  is  only  dependent  on  DPF  system  parameters. 

Since  F/W  ratios  are  still  below  the  desired  value,  it  is  obvious  that  extra  thrust  is  needed 
while  keeping  system  mass  down.  This  may  be  possible  by  adding  additional  thrusters  if  the 
extra  thrust  produced  more  than  compensates  for  the  extra  mass  due  to  capacitors,  shielding  etc. 
The  additional  thrusters  will  not  affect  1^  values  and  Figure  21  is  still  valid.  Adding  additional 
thrusters  is  also  insurance  against  a  failure  in  any  one  thruster.  Figures  22  and  23  show  the 
general  trends  of  vehicle  F/W  as  additional  thrusters  are  added  for  10  and  40  km/s  mission  Av’s 
respectively.  For  low  Av  missions,  where  required  propellant  mass  is  small,  it  is  best  to  run 
several  thrusters  at  high  propellant  mass  flow  rates.  Figure  23  best  to  run  several  thrusters  at 
relatively  low  propellant  mass  flow  rates.  Initial  mass  in  low  earth  orbit  (IMLEO)  is  probably 
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Figure  15.  Av  Capability  vs.  Payload  Mass  Fraction 
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Figure  16.  Vehicle  F/W  vs.  Propellant  Mass  Flow  Rate  for  Av=5  ktn/s 
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Figure  19.  Vehicle  F/W  vs.  Propellant  Mass  Flow  Rate  for  Avs30  km/s 
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Figure  20.  Vehicle  F/W  vs.  Propellant  Mass  Flow  Rate  for  Av=40  km/s 
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Figure  32.  F/W  Ratio  vs.  Propellant  Mass  Flow  Rate  and  Number  of  Thrusters  for  Av=10  km/s 
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the  limiting  factor  fm*  the  number  of  thrusters  since  each  thruster  increases  system  mass  and  it  is 
quite  expensive  to  raise  payload  into  LEO.  Figure  24  shows  quantitative  results  for  20  MA  and 
10. 20, 30  and  40  km/s  Av’s,  respectively.  Figure  23  shows  that  F/W  ratios  approaching  0.075  at 
an  Itp  of  about  4000  s  are  possible  with  4  thrusters  and  a  propellant  mass  flow  rate  of  4  kg/s. 

This  value  seems  to  be  approaching  those  necessary  for  a  manned  Mars  mission.  In  this  mode, 
the  DPF  propulsion  system  can  be  competitive  with  chemical  and  nuclear  fission  rockets  for  use 
in  distant  space  missions. 

It  is  desirable  to  have  F/W  ratios  as  large  as  those  offered  by  other  means  of  propulsion, 
namely,  chemical  and  nuclear  fission  rockets.  At  this  time,  these  F/W  values  seem  to  be  largo* 
than  the  DPF  propulsion  system  can  produce.  However,  engineering  breakthroughs  and 
innovations  may  make  these  levels  attainable.  If  capacitor  specific  energies  can  be  increased  to 
about  20  kJ/kg,  system  masses  may  be  reduced  until  high  F/W  ratios  are  possible.  Improvement 
in  confinement  time  or  repetition  rate  would  also  serve  to  increase  F/W  by  increasing  the  exhaust 
power. 

These  parametric  calculations  are  very  sensitive  to  the  assumed  dimensions  of  the  pinch. 
Since  the  fusion  power  depends  on  the  plasma  volume,  any  error  in  the  estimation  of  pinch 
dimensions  may  greatly  understate  the  amount  of  fusion  power  produced.  This  will  in  turn 
affect  both  thrust  and  Lp.  Caution  should  be  taken  when  "using"  these  numbers  in  the 
realization  that  final  results  can  depend  greatly  on  assumed  initial  parameters.  A  computer  code 
which  calculates  important  propulsion  parameters  given  the  initial  DPF  parameters  (e.g.. 

Table  3)  can  be  found  in  Appendix  B. 


II.2.  PARAMETRIC  STUDIES  OF  DPF  WITH  SPIN-POLARIZED  D-3He  FUEL* 

Recently,  the  effect  of  nuclear  spin  polarization  of  /usion  fuels  on  fusion  reactor  operations 
has  been  investigated  (48).  The  plasma  focus  device  a  nalyzed  Section  II.l  is  based  on  the 
"Livermore-l"  dense  plasma  focus  without  spin-polarization.  Table  4  lists  all  the  assumed 
parameters  used  in  the  evaluation  including  the  capacitor  banks  being  capable  of  delivering  a 
maximum  current  in  excess  of  10MA. 


*  The  major  portion  of  this  section  of  the  report  has  been  A-cumented  as  the  M.S.  project  report  by  M. 
Wang  as  a  partial  requirement  for  the  MS.  degree  from  Punk*  University. 
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Figure  24.  Vehicle  F/W  vs.  Propellanl  Mass  Flow  Rate  for  1*20  MAantLA  v  =  10  40  kvrt$ 
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Consistent  with  the  earlier  analysis  in  Section  HI,  the  impulsive  operation  mode  is 
considered  to  be  roost  suitable  for  interplanetary  travel  due  to  hs  higher  F/W  ratios.  Thus  this 
section  attempts  to  make  a  comparison  for  impulsive  operation  of  DPF  propulsion  system  with 
the  spin-polarized  fuel.  Since  the  system  mass  is  vexy  important  for  fusion  propulsion  with  the 
capacitor  mass  being  the  major  one  among  the  system  components,  the  scaling  of  die  capacitor 
mass  with  its  current  becomes  a  key  issue  in  die  analysis.  For  the  "LivermoreT"  device,  it  was 
initially  assumed  that  the  capacitor  mass  scaled  as  the  current  squared.  However,  as  indicated 
earlier  in  Section  II.  1  and  also  in  Appendix  A,  it  can  be  shown  that  the  capacitor  mass  might 
possibly  scale  as  I®3.  In  this  section,  both  cases  are  discussed.  Multiple  thrusters  are  also 
considered  since  they  not  only  provide  insurance  against  a  failure  but  also  can  raise  the  F/W 
values.  The  extra  thrust  produced  by  additional  thrusters  more  than  compensates  for  die  extra 
mass  due  to  capacitors,  shielding,  etc.  As  shown  earlier  in  Figs.  22  and  23  for  I2  scaling,  die 
F/W  values  with  scaling  (e.g..  Figures  42  and  43)  are  more  pessimistic  than  those  with  I2 

scaling,  and  in  this  case,  multiple  thrusters  are  much  less  effective  than  for  I2  case. 

Figure  26  for  the  I2  scaling  case  shows  the  variation  of  specific  impulse  for  any  Av  as  well  as 
any  thruster  number  since  the  propellant  exit  velocity  is  only  dependent  on  DPF  system 
parameters  and  the  additional  thrusters  will  not  affect  Ijp  values.  From  this  figure,  it  is  obvious 
that  the  spin  polarization  can  increase  the  specific  impulse.  Figures  27  and  28  illustrate  the 
dependence  of  F/W  on  propellant  mass  flow  rate  and  F/W  vs.  1^,,  respectively,  for  low 
Av  (Av  =  10  km/s)  from  1  to  4  thrusters.  Figures  29  and  30  show  the  same  parameters  for  high 
Av  (Av  =  40  km/s).  All  the  figures  prove  again  that  the  spin  polarization  makes  the  F/W  values 
increase.  Figure  27  seems  to  show  that  F/W  can  be  increased  indefinitely  by  simply  exhausting 
more  propellant;  however,  this  lowers  the  propellant  outlet  temperature  which  in  turn  increases 
the  total  amount  of  propellant  required  for  the  mission.  This  is  not  apparent  in  missions  with 
small  Av’s  since  the  originally  required  propellant  mass  is  so  small  compared  to  other  system 
masses.  Figure  29  actually  shows  the  cases  where  propellant  becomes  the  dominant  mass  and  an 
increase  in  propellant  mass  flow  actually  decreases  system  F/W.  But  as  one  can  see  from 
Figures  27  and  29,  the  highest  current  no  longer  leads  to  the  highest  F/W  as  number  of  thruster 
increased.  The  optimum  current  now  occurs  at  about  20MA.  This  is  because  for  more  thrusters 
the  more  propellant  mass  must  be  carried  along  with  capacitor  mass  at  higher  current.  And  since 
the  spin  polarization  can  reduce  the  propellant  mass,  this  makes  the  shift  of  the  curves  at  20MA 
and  25MA  more  obvious.  The  crossover  between  the  20  and  25MA  case  is  due  to  the 
additional  capacitor  mass  being  compensated  for  by  the  additional  propellant  mass  which  must 
be  carried  along  at  the  lower  current.  The  x-axes  of  Figures  28  and  30  indicate  propellant  mass 
flow  decreasing  while  specific  impulse  increases  according  to  specific  impulse  changing 


47 


CIFIC  IMPULSE  <SEC ) 


16500 

14850 

13200 

11550 

9900 

8250 

6600 

q-  4950 

t n 

3300 

1650 

0 

0 


WLH1CUE  JiWlJSl-TO-UCIGMI  RATIO 


■  ifM  (•mhmIvim) 

•  »an» 

•  2tn* 


(*in-Mlir  iM) 


17  30  99  <1  V*  C3  OO 

P9:?ELLANT  hASS  rL0U  «OT£  (KG-'S) 

43  Ba  *  II 

moss  rLSU  ROT! 

a)  1  thruster 

b)  2  thrusters 

10"*  (Min*t«l«r  ii*4> 

1 V*  (MinftliritM) 

W*»  !W)n*HlrilN) 
J9**  I  Win*Ml«r  ntci 
)#">  (tmiir.tN) 


go  jooj  j-  W 


10  is  oo  eo  oo  M  oo 

propeii.ont  moss  tlcw  ro*:  (kOS> 


5  00  [0  30  19  00  JO  30  J9  30 

PR0PEO.ANT  MOSS  "L3U  ROTE  ;kG'S> 


c)  3  thrusters 


d)  4  thrusters 


Figure  27.  Vehicle  F/W  vs.  Propellant  Mass  Flow  Rate  for  Av=10km/s  (I2  Scaling) 
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Figure  30.  Vehicle  F/W  vs.  Specific  Impulse  for  Av=40km/s  (I2  Scaling) 
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adversely  with  propellant  mass  flow  as  shown  in  Figure  26.  The  curves  bend  differently  because 
the  F/W  increases  with  propellant  at  low  Av  but  has  an  optimum  value  at  propellant  mass  flow 
about  5  kg/s  (from  Figure  29),  i.e.,  specific  impulse  about  4000  s  (from  Figure  26)  for  spin- 
«  polorized  case.  Since  the  optimum  current  is  20MA  and  F/W  values  increase  with  number  of 

thrusters,  Figure  31  and  32  demonstrate  the  trend  of  F/W  vs.  propellant  mass  flow  and 
F/W  vs.  Isp  for  different  Av,  Obviously,  DPF  propulsion  can  produce  higher  F/W  for  low  Av 
missions.  This  is  consistent  with  the  implication  in  Sect  n.l  that  DPF  propulsion  systems  are 
best  suited  for  low  Av  missions. 

For  the  I8/3  scaling  case,  almost  the  exact  same  information  is  presented  except  that  the  F/W 
values  are  lower  and  the  optimum  current  becomes  15MA  as  shown  in  Figures  33  to  36.  This  is 
due  solely  to  the  increasing  capacitor  mass  at  higher  currents  which  was  not  so  evident  before 
because  of  the  I2  scaling.  Comparing  Figure  24  to  Figure  33,  and  Figure  35  to  Figure  36,  one 
notes  thruster  performance  being  adversely  affected  by  the  higher  capacitor  masses  due  to  the 
new  lm  scaling  relation.  The  increasing  number  of  thrusters  does  little  to  increase  F/W  and  the 
raising  of  the  propellant  mass  flow  rate  drastically  reduces  F/W  because  of  the  large  propellant 
masses  involved.  Again,  spin  polarization  reduces  propellant  mass  requirements  at  higher  Av 
which  makes  the  15MA  curve  rise  up  more  rapidly. 

From  the  analysis  above,  one  recognizes  that  there  is  a  corresponding  optimum  condition  for 
each  Av  and  scaling  relation.  The  conditions  which  can  provide  the  maximum  F/W,  i.e.,  Av  *  10 
km/s,  I  ■  20MA,  and  using  4  thrusters,  are  chosen  to  make  the  comparison  between  the  base  case 
and  the  spin-polarized  case  as  listed  in  Tables  5-1  and  5-2  for  low  and  high  Av’s,  respectively. 
Various  parameters  in  Tables  5-1  and  5-2  are  described  in  the  subsequent  discussions. 

The  total  fusion  power  produced  by  DPF  is 

PF  =  K  tplp [nDnHe  <  OV  >  DH,  WDH<  +  i-n&  («JV>DDnWDDt  +  <0V>DDp  WDDf)]Rrept  (17) 

,  where  Rrep  is  number  of  firings  per  unit  time  (s-1 ),  and  t  is  pinch  stable  time(s),  as  described  in 

Table  4.  Due  to  low  cross-sections  at  the  operating  temperatures  of  D-3He,  secondary  reactions 
including  D-T ,  T-T,  and  3 He-3 He  reactions  are  not  included  in  Eq.(17). 


53 


Aw  »  tOKivg 
Aw  •  30  KM/’S 
Aw  •  30KIVS 
Aw  •  HOKM/S 
Aw  ■  10KIUS 
Aw  »  30KM/S 
Aw  *  30KM/S 
Aw  ■  WOKtV$ 


<spirv>pol«r’iz»d) 

<  >P  i  rvpo  l  «p  i  z  tg  ) 

<*pin-ool«riz*p) 

<*Pln-pol«riz»<3) 

<wnpol«r iz*d) 
<unpol4r  Izad) 
<un0ol«r  tz»<j> 
<unool*r ix*d) 


_L 


10.00 


13.00 


propellant "mass  fLO? RATS 


25.00 


30. 00 


vi  "ss®  for  av = 10  «  «*  **■»*> 


urn i n  r  thrust- to -ucight  ratio 


0,7000 


0  *6000 


,9000 


0  .4000 


0 ,3000 


0 .2CCC 


0  ,1000 


0  .0000 


Aw  •  10KIVS  C  spin-polar iztd  > 
Aw  «  20KH/S  <  spin-polar iztd  > 
Aw  »  30KIV$  <  spin-polar iztd  > 
Aw  *  40KM/S  <  spin-polar iztd  > 
Aw  »  lOKflrS  <  unpo l ar iztd  > 

Aw  «  20KMr$  <  unpolariztd  > 

Aw  •  30Kf1/S  <  unpolariztd  > 

Aw  *  40KP1/S  <  unpolariztd  ) 


Q.  2000  .  4000.  6000,  8000,  10000  ,  12C00  .  14000 

SPECIFIC  IMPULSE  <SEC> 


Figure  32.  Vehicle  F/W  vs.  Specific  Impulse  for  Av  (I2  Scaling) 
(I=20MA,  4  thrusters) 


55 


MiMV  IMHO  MMIfU  IHHUSI  10  urtttll  «ftllO 


Table  5.1.  Propulsion  Parameters  for  Base  Cast  afcd  Spirt-Pdiarfred  Case 
( Av  =  lOkrn/s,  1  =  20  MAy  *'  =  30kg/s) 


Base  case  ,  ....  Spifcl^&ifrgd_„ 


D  Bumup  Fraction 

fo 

0.63 

0.70 

0.63 

070 

3  He  Bumup  Fraction 

0.48 

0.58 

0.48 

OSfe 

D-3He  Fusion  Power 

Pbtk 

(MW) 

2959.41 

3582.31 

2939.41 

358231 

DD„  Fusion  Power 

PDDn 

(MW) 

11.09 

9.05 

11.09 

9.05 

DDp  Fusion  Power 

Poop 

(MW) 

49.66 

40.50 

49.66 

4&30 

Total  Fusion  Power 

Pp 

(MW) 

3020.15 

3631.86 

3020.15 

3631.86 

Power  to  Focus 

Pm 

(MW) 

88.10 

88.10 

223.43 

223.43 

Bremsstrahlung  Loss 

Pb 

(MW) 

16.28 

13.09 

16.28 

13.09 

Cyclotron  Loss 

Pc 

(MW) 

84.23 

74.42 

84.23 

74.42 

Total  Power  Lou 

Pl 

(MW) 

100.51 

87.51 

100.51 

87.51 

Power  Increase 

AP 

(MW) 

2831.54 

3456.25 

2696.20 

3320.92 

Total  Mass  Flow 

Mr 

(kg/s) 

31.02 

30.96 

31.02 

30.96 

Propellant  Thrust 

fp 

(N) 

4.85X105 

5.25x10s 

4.85x10s 

5.25x10s 

Total  Bum  Time 

lb 

(s) 

4852.44 

4360.00 

1.56x10* 

1.40x10* 

Payload  Mass 

Ml 

(kg) 

lxlO5 

1x10s 

1x10s 

1x10s 

Propellant  Mass 

Mp 

(kg) 

1.48x10s 

1.33x10s 

4.76x10s 

4.27x10s 

Prop.  Sys,  Mass 

Mpiyt 

(kg) 

2.22x1 04 

1.99x10* 

7.14x10* 

6.40x10* 

Fuel  Mass 

Mp 

(kg) 

1.96 

1.76 

6.29 

5.65 

Fuel  Sys.  Mass 

Mpm 

(kg) 

0.20 

6887.88 

0.63 

4.43x10s 

Capacitor  Mass 

W 

Me 

(kg) 

6.89x10* 

6.89x10* 

4.43x10s 

4.43x10s 

Shield  Mau 

M* 

(kg) 

1.27x10* 

1.30x10* 

1.37x10* 

1.40x10* 

Magnet  Mass 

Mb 

(kg) 

270.20 

270.20 

270.20 

270.20 

Total  Mass 

Mr 

(kg) 

3.52x10s 

3.42x10s 

1.10x10s 

1.09x10s 

Total  Thrust 

F 

(N) 

1.94x10s 

2.10x10s 

1.94x10s 

2.10x10s 

Thrust-to-Weight 

F/W 

0.56 

0.63 

0.18 

0.20 

Specific  Impulse 

w 

(s) 

1622.50 

1760.84 

1622.50 

1760.84 

Increased  Factor  of  Fusion  Power 

1.20 

1.20 

Increased  Factor  of  AP 

1.22 

1.23 

Increased  Factor  of  F/W 

1.13 

1.11 

Increased  Factor  of  1» 

1.09 

1.09 

*  Additional  hydrogen  used  for  propellant  besides  coolant  mass  flow. 
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Table  5.2.  IVotwdriow  PnrMMtCT  for  Baae  Case  and  Spin  WunedCwe 

(Av  =  4Okmfs,I  =  20MA,M-*  =  i*kg/s) 


Base  case  Spin-Polarized 


D  Bumup  Fraction 

0.63 

0.70 

0.63 

0.70 

3  He  Bumup  Fraction 

fife 

0.48 

0.58 

0.48 

0.58 

D-3He  Fusion  Power 

PDHe 

(MW) 

2959.41 

358131 

2959.41 

3582.31 

DDn  Fusioi.  Power 

PDDn 

(MW) 

11.09 

9.05 

11.09 

9.05 

DDp  Fusion  Power 

PDDp 

(MW) 

49.66 

40.50 

49.66 

40.50 

Total  Fusion  Power 

Pf 

(MW) 

3020.15 

3631.86 

3020.15 

3631.86 

Power  to  Focus 

Pin 

(MW) 

88.10 

88.10 

223.43 

223.43 

Bremsstrahlung  Loss 

Pb 

(MW) 

16.28 

13.09 

16.28 

13J09 

Cyclotron  I  oss 

Pc 

(MW) 

84.23 

74.42 

84.23 

74.42 

Total  Power  Loss 

Pl 

(MW) 

100.51 

87.51 

100.51 

87.51 

Power  Increase 

AP 

(MW) 

2831.54 

3456.25 

2696.20 

3320.92 

Total  Mass  Flow 

Mr 

(kp's) 

5.02 

4.96 

5.02 

4.96 

Propellant  Thrust 

Fp 

«N) 

1.72x10s 

1.88x10s 

1.72x10s 

1.88x10s 

Total  Bum  Time 

tb 

(s) 

6.95x10* 

6.01x10* 

2.23x10s 

1.93x10s 

Payload  Mass 

ml 

(kg) 

1x10s 

1x10s 

1x10s 

1x10s 

Propellant  Mass 

Mp 

(kg) 

3.14x10s 

2.68x10s 

l.OlxlO6 

8.60x10s 

Prop.  Sys.  Mass 

Mp.y, 

(kg) 

4.71x10* 

4.02x10* 

1.51x10s 

1.29x10* 

Fuel  Mass 

mf 

(kg) 

28.02 

24.22 

89.99 

77.80 

Fuel  Sys  Mass 

Mpjy, 

(kg) 

6890.13 

9.00 

4.43x10* 

Capacitor  Mass 

Me 

(kg) 

6.89x  i  1 

6.89x10* 

4.43x10s 

4.43x10s 

Shield  Mass 

Msh 

(kg) 

1.51x10 

’.53x10* 

1.62x10* 

1.64x10* 

Magnet  Mass 

7lB 

(kg) 

270.20 

270.20 

270.20 

270.20 

Total  Mass 

Mt 

(kg) 

5.45x10s 

4.99x10s 

1.72X106 

1.59xl06 

Total  Thrus* 

F 

(N) 

6.89X106 

7.52x10s 

6.89x10s 

7.52x10s 

Thrust-to-Weight 

F/W 

0.13 

0.15 

4.09xl0"2  4.81xl0-2 

Specific  Impulse 

I*p 

(s) 

3888.86 

4300.36 

3888.86 

4300.36 

Increased  Factor  of  Fusion  Power 

1.20 

1.20 

Increased  Factor  of  AP 

1.22 

1.23 

Increased  Factor  of  F/W 

1.15 

1.18 

Increased  Factor  of  1^ 

1.11 

1.11 

♦Additional  hydrogen  used  for  propellant  besides  coolant  mass  flow. 
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Hie  initial  number  density  1*0=0^  =  ynj  =  1.18  x  1026  m-3  and  die  reaction  rate  parameters 

of  <ov>^|_.  <ov>nn_.  and  <qv>nn^  arc  used,  consistent  with  the  values  used  for  the 
conventional  D-3He  case  in  Section  IL1.  For  spin-polarized  case,  «tv>d^  =  l^<ov>Dkfc> 
<ov>poii  =<ov>m)>.,  and  <ov>|)|^  =<ov>ix]^  were  used.  During  the  pinch  stable  time,  the 
number  densities  of  D  and  3 He  were  actually  functions  of  time  which  decreased  with  time. 

Since  die  D-3He  reaction  rate  was  increased  by  spin  polarization  and  the  deuterium  panicles 
decreased  at  faster  rate,  the  fusion  power  produced  by  DDn  and  DDp  reactions  were  slightly 
decreased.  Deuterium  bum  up  fraction  was  0.63  as  compared  to  that  of  3He  with  0.48  as  listed 
in  Table  5. 

As  one  can  see  from  Tables  5.1  and  5.2,  the  total  fusion  power  was  increased  by  a  factor 
about  1.2.  Also  one  can  find  that  brcmsstrahlung  and  synchrotron  losses  are  slightly  suppressed 
by  spin  polarization  because  ot  the  faster  decrease  of  D  and  3He  number  densities.  The 
increased  fusion  power  and  the  decreased  total  power  loss  lead  to  the  fact  that  more  power  can 
be  absorbed  by  propellant  to  produce  the  thrust  for  spin-polarized  case.  This  is  why  the  total 
thrust,  the  sum  of  prc^cSant  thrust  and  the  thrust  produced  by  expelled  fuel  increased.  One 
thing  should  be  noted  is  that  the  power  required  to  operate  plasma  focus  is 

Pin=IVRreptdis,  (18) 

where  t^j  is  time  for  fill  gas  to  be  discharged  with  other  parameters  in  Eq.(18)  given  in  Table  4. 
For  the  I2  scaling  case,  V  «  i,  then  « I2.  But  for  the  I8/3  scaling  case,  V  «  I4/3,  then 
Pi„  oc  l7/3.  Thus  this  current  relationship  is  very  sensitive  to  the  different  scaling  laws. 

Part  of  the  synchrotron  radiation  will  be  absorbed  in  walls  and  electrodes.  Since  the 
synchros  was  lowered  by  spin  polarization,  the  amount  of  power  that  should  be  removed  from 
the  walls  and  electrodes  is  correspondingly  decreased.  Therefore  the  coolant  mass  flow  which  is 
needed  to  keep  inlet  temperature  to  turbine  below  2,000°K  is  lowered;  hence,  the  total  mass 
flow  which  is  the  sum  of  additional  propellant  and  coolant  mass  flows  is  reduced  for  the  spin- 
polarized  case. 

Since  the  specific  impulse  is  increased  due  to  increased  fusion  power,  the  exhaust  velocity, 
vex,  is  also  increased.  Thus  the  propellant  system  masses  are  decreased  via 

Mp  =  (eAvAl~  -  1)  (Ml  +  Mc  +  MB) ,  (19) 
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assuming 


Mpsy,  =0.15Mp . 


(20) 


Reduced  propellant  mass  also  reduces  the  propellant  bum  time,  and  thus  the  fuel  and  fuel  system 
masses  are  reduced  since  Mp  « tb  and  =  O.IMp. 


As  shown  in  Tables  5.1  and  5.2,  capacitor  mass  for  scaling  case  is  much  bigger  than  that 
for  I2  scaling  since  at  I8/3  scaling 


Me 

Mco 


r  > 

V 

2 

r  > 

I 

V0 

Io 

to  to 

(21) 


and  this  large  mass  makes  the  F/W  ratios  in  I8/3  scaling  much  smaller  than  those  of  I2  scaling. 

The  thrust-to-weight  ratios  in  Tables  5.1  and  5.2  increase  due  to  the  raised  thrust  as  well  as 
lowered  mass,  and  the  specific  impulse  also  increased  due  to  increased  fusion  power  by  the  spin 
polarization  for  both  I2  and  18/3  scaling  cases. 
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m.  CONCLUSIONS  AND  RECOMMENDATIONS 


If  plasma  temperature  scaling  holds,  the  dense  plasma  focus  could  be  a  relatively  easy  way 
of  obtaining  hot,  high  density  plasmas.  Further  study  of  the  dynamics  of  the  pinching  process 
and  of  the  eventual  disruption  of  the  pinch  due  to  MHD  instabilities  is  necessary.  Trapping  a 
strong  axial  magnetic  field  inside  the  pinch  may  stabilize  it  by  reducing  the  rate  at  which  the 
instabilities  grow.  Accurate  measurements  of  pinch  dimensions  is  critical  in  the  computation  of 
important  propulsion  parameters. 

The  most  important  issue  to  be  resolved  is  how  various  parameters  scale.  Scaling  of  plasma 
temperature  with  capacitor  current  must  be  investigated  at  currents  larger  than  1  MA,  where  I2 
scaling  may  break  down.  Also  of  extreme  importance  is  the  scaling  of  capacitor  mass  with 
current.  Both  I2  and  I8/3  scalings  were  considered  in  this  study. 

For  impulsive  thrusting  for  non-spin  polarized  case,  the  increase  in  capacitor  mass  is 
sufficiently  large  to  reduce  F/W  values  beyond  currents  of  20  MA  for  impulsive  thrusting.  The 
scaling  also  determines  the  effectiveness  of  using  multiple  thrusters  on  a  mission.  At  I2  scaling, 
multiple  thrusters  serve  to  greatly  increase  F/W  values  as  well  as  serve  as  insurance  against  the 
failure  of  a  thruster.  However,  at  I8(3  scaling,  multiple  thrusters  are  much  less  effective  in 
increasing  F/W  because  of  the  large  increase  in  capacitor  mass.  For  Av  =40  km/s,  F/W  ratios  of 
almost  0.08  are  possible  at  around  4000  sec  of  Isp  with  4  thrusters  at  20  MA.  This  seems  to  be 
the  optimum  operating  regime  for  the  DPF  propulsion  system  in  the  impulsive  firing  mode. 

The  F/W  ratios  for  the  continuous  firing  modes  of  operation  of  the  DPF  propulsion  system 
for  non-spin  polarized  case  fall  about  2  orders  of  magnitude  short  of  that  which  is  required  for 
manned  space  travel.  Although  none  of  the  three  modes  of  operation  of  the  DPF  propulsion 
system  are  suitable  for  long  missions  in  space,  each  may  have  a  mission  for  which  it  is  very  well 
suited.  They  are  probably  best  suited  for  low  Av  missions  such  is  orbital  transfer  or  perhaps  a 
lunar  shuttle,  which  do  not  require  large  F/W  values. 

The  greatest  benefit  from  spin  polarizing  D-3He  analysis  is  an  increase  in  the  D-3He  output 
power  by  20%  up  to  possibly  50%.  The  power  increase  by  spin  polarization  indeed  increases  the 
thrust-to- weight  ratio  and  specific  impulse.  With  the  possible  suppression  of  D-D  reactions,  this 
50%  power  increase  could  actually  be  envisioned  as  the  total  power  increase.  This  indicates  that 
a  50%  power  increase  with  little  or  no  reactor  modifications.  Examining  the  benefit  to  fusion 
propulsion  one  finds  the  actual  total  power  increase  in  less  than  50%  because  of  the  decrease  of 
ion  number  densities  during  the  pinch  lifetime,  which  leads  the  reaction  rate  to  decrease  with 
time.  The  spin-polarized  case  presented  here  can  probably  be  improved  by  optimization  of  the 
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plasma  and  the  propulsion  parameters;  e.g.,  this  work  treated  the  DD  reaction  as  unchanged  by 
spin  polarization.  Though  the  F/W  values  and  specific  impulse  indeed  increased  by  adopting 
spin  polarization,  for  high  Av  missions,  however,  F/W  ratios  of  the  DPF  propulsion  system  is 
still  too  low  (cf.  0.2).  At  present,  DPF  propulsion  is  best  suited  for  low  Av  missions  as  described 
above. 
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APPENDIX  A 


Is0  CURRENT  SCALING  RELATION 


It  is  important  to  get  good  estimates  of  system  masses  in  older  to  make  accurate  calculations 
of  system  F/W  ratios.  One  very  important  mass  to  be  calculated  is  the  capacitor  mass.  Itis 
necessary  to  extrapolate  the  required  electrical  parameters  from  some  reference  case.  The 
reference  case  is  chosen  as  the  Uvennoie-r  device  and  is  denoted  here  by  the  subscript  "o". 
Assuming  a  constant  capacitance  and  capacitor  specific  energy,  die  capacitor  charging  potential 
can  be  found  as  a  fu.iction  of  current  by  consistently  solving  Eqs.  7-10  with  Eq.  2  in  Section 
ILL  It  can  be  shown  that  the  charging  potential  scales  rs  I4** 

V_ 

V0 


Since  capacitor  mass  is  proportional  to  capacitor  energy  which  is  proportional  to  the  square  of 
the  charging  potential,  capacitor  mass  scales  as  I8/3 


Me 

Mc,o 


Figures  37-4 .  show  F/W  ratios  for  various  Av’*  for  the  new  scaling  law.  Immediately 
apparent  in  these  Figures  is  that  the  highest  current  no  longer  leads  to  the  highest  F/W.  This  is 
due  solely  to  the  increasing  capacitor  mass  at  higher  currents  which  was  not  evident  before 
because  of  the  I2  scaling.  The  optimum  current  now  occurs  at  about  20  MA.  Figures  39  and  40 
show  a  crossover  between  the  IS  and  25  MA  cases.  It  is  here  where  the  additional  capacitor 
mass  is  compensated  for  by  the  additional  propellam  mass  which  must  be  carried  along  at  the 
lower  current.  Figures  39-41  show  that  F/W  decreases  as  additional  propellant  flow  is  added  due 
to  the  fact  that  propellant  is  now  the  dominant  system  mass.  It  should  be  recognized  that  at 
some  very  high  propellant  mass  flow  rate  (or  at  a  Jgher  Av)  the  highest  current  will  eventually 
yield  the  highest  F/W  ratios  as  the  difference  in  capacitor  mass  is  small  compared  to  propellant 
mass. 


Again,  it  is  possible  to  use  multiple  thrusters  to  enhance  performance.  Figures  42  and  43 
show  thruster  performance  at  20  MA  and  10  and  40  km/s,  respectively.  Thruster  performance  is 
obviously  adversely  affected  by  the  higher  capacitor  masses  due  to  the  new  scaling  law 
(compare  to  Fig.  22),  and  the  additional  thrusters  have  only  a  minor  effect  on  vehicle  F/W  for 
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low  Av  while  significantly  increasing  IMLEO.  Similarly,  in  Figure  43,  the  effect  of  die 
increasing  number  of  thrusters  on  F/W  is  not  as  significant  as  for  die  I2  scaling  (e.g.,  Figure  23). 
Raising  the  propellant  mass  flow  rate  drastically  reduces  F/W  because  of  thefcuge  propellant 
masses  involved.  Figures  44  and  45  give  some  quantitative  results  and  shows  (hat  for  Av  =  40 
km/s  the  DPF  propulsion  system  is  still  about  in  oider  of  magnitude  below  desirableT/W  ratios 
for  a  manned  Mars  mission. 


76 


uni  1 1  cle  ti  irijst-to  -he  ir.i it 


0.00  2  00  10  00  IS. 00  20.00  25.00  3C  .00  35 

FRO PELLAINT  nASS  FLOW  RATE  (KG'S) 


.  Figure  37.  Vehicle  FAV  vs.  Propellant  Mass  Flow  Rate  for  A  v  =  5  km/s  (I8/3  scaling) 
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Figure  38.  Vehicle  F/W  vs.  Propellant  Mass  Flow  Rate  for  A  v  =  10  km/'s  (I8/3  scaling) 
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Figure  40.  Vehicle  F/W  vs.  Propellant  Mass  Flow  Rate  for  A  v  =  30  km /s  (I8/3  scaling) 
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Figure  41.  Vehicle  F/W  vs.  Propellant  Mass  Flow  Rate  for  Av  =  40km/s  (18/J  Scaling) 
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Figure  42.  F/W  Ratio  vs.  Propellant  Mass  Flow  Rate  and  Number  of  Thrusters 
for  Av  =  lOkm/s  (r3  Scaling) 
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Figure  44.  Vehicle  FAV  vs.  Propellant  Mass  Flow  Rate  for  1=20  MA  and 
Av  =  lOkm/s  (I8/3  Scaling) 
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APPENDIX  B 


CODE  FOR  EVALUATION  OF  DPF  PROPULSION  WITH 
NON-POLARIZED  AND  SPIN-POLARIZED  FUELS 

How  to  use  the  code: 

This  is  a  simple  direction  for  how  to  use  the  DPF  fusion  propulsion  code  with  D~3He  fuel. 

1.  There  are  four  major  programs,  all  of  them  are  for  impulsive  operation: 

a)  impsq.f :  I2  scaling,  unpolarized  fuel 

b)  spimpsq.f :  l2  scaling,  spin-polarized  fuel 

c)  imp83.f :  1 8/3  scaling,  unpolarized  fuel 

d)  spimp83.f :  Iw  scaling;  spin-polarized  fuel 

2.  Input  file :  test.dat 

3.  Compile  the  program  as  usual,  i.e.,  T77  impsq.f  or  'fort  impsq.f 

4.  Use  ‘a.out’  to  get  the  output  files.  Each  program  produces  four  output  files: 

a)  Test.out :  a  file  lists  all  the  information  you  need 

b)  tw.out :  used  to  plot  thrust-to-weight  ratio  vs.  propellant  mass  flow  rate  graph  with 
genplot 

c)  isp.out :  I,p  vs.  propellant  mass  flow  rate 

d)  twisp.out :  thrust-to-weight  ratio  vs.  lv 

5.  The  program  imp3d.f  is  just  for  reference.  It  can  be  used  to  create  3-D  plot  with  PLOT3D. 
One  thing  should  be  noticed:  modify  it  before  you  use  to  avoid  any  possible  error. 

6.  Whenever  you  rerun  the  program,  remember  to  remove  all  four  output  files  first. 


nnnnn  oon  oooo 


C 

Q  *★***#**★  *#•#★★* Hr* ************ *#*w*#*  Hr* w***#**r  -****»**#i**4*W*f*****’*,#***A 

C  *********  THE  DENSE  PLASMA  FOCOS  :  A  FUSION  PROPULSION  SCENARIO  &/***&?.  ,, 

£  *********  *t*#*w^tt*f* ********** v*******#********************#**#*****^#*^ 

*  **************** ************************* ****** * ******* 

#  i  .  /  ••  Jf  i  4 


PROGRAMMER  :  CHRISTOPHER  L.  LEAKEAS 

SCHOOL  Or  NUCLEAR  ENGINEERING 

PURDUE  UNIVERSITY 

WEST  LAFAYETTE,  IN  47907 

WRITTEN  AT  :  ASTRONAUTICS  LABORATORY 
AL (AFSC) /LSVF 

EDWARDS  AFB,  CA  93523-5000 
MODIFIER  :  ME I -YU  WANG 

THIS  PROGRAM  WILL  CALCULATE  THE  PARAMETERS  NECESSARY  IN  THE  OPERATION 
OF  A  DENSE  PLASMA  FOCUS  FOR  USE  AS  A  SPACE  THRUSTER. 


PROGRAM  FOCUS 

RIAL  ADMSFLW,  AVMASS , BPNCK, BCLTZ, CPH2, CPELEC, CPHION, CAP 
REAL  CONST, CYCREFL, DFRACT, DSCKRG. DHE1 , DHS2.DHE3 , DHE4 
REAL  DDN1 , DDN2 , DDN3 . DDN4 , DDP1.DDP2, DDP3 , DDP 4 
REAL  ENERGY, ELECTHR, ELECTEN, F, FRACT, FSNPLW, FPNCH 
REAL  GRAV,  HIONTHR,  HEFRACT,  IMAXSQ,  IMAX,  ISP,  IVCL 
REAL  IMAGNET, IMOPT, KT , KTOPT , LANODE ,  LPNCH,  LIMIT 
REAL  MC AP ,  MP ROP ,  MPROSTR,  MMAGNET ,  MP K'i LD ,  MFUEL ,  MFUELS Y 
REAL  MSHIELD ,  MSFLW,  MTOT, MUNOT, MH2 , NFTRRST, PNCHRAD ,  P VOL 
REAL  PNCHTIM, PABSW, PI2R, PIN, PI, PMAGNET,  PROPTHR, QDOTREM 
REAL  QLEFT, RHOI, RA, RC, RHOCU, REPRATE, SIGVDHE 
REAL  SIGVDDN, SIGVDDP, SPECEN, TBURN, TOTHRST, TF? , TPLOSS 
REAL  TEMP . TDSLTAP , TPFDHE, TPFDDN, TPFDDP , TAVG, TP3REM, TPCYC 
REAL  THICK, TDI5ION, TSTAGH, TSTAGE, TTHRTH,  TTHRTE, VOLT, VRUN 
REAL  7KI0NSX, VHICNTH, VELECTH, VELSCEX, WEHS, WCDN, WDDP, WASTE 
REAL  X, X3ECTAR, VEX. MIN, VMAX, KNOT 

C  DECLARE  ARRAYS  FOR  ITERATIONS 

REAL  DELTA?  {1000)  ,DNP(10Cl)  ,  KEN?  (1001) ,  NPNCKC0C1) 

REAL  PFDKE ( 1000 ) , PFDDN ( 1000 ) . PFDDP ( 1000 )  ,  PFTC7 ( 10C0 ) 

REAL  P9REM (1000) .PCYC(lOOO) , PLOSS (100C ) , RRDKZ (1000 ) 

REAL  RRDDN(ICOO)  .RRDDP(IOCO) 

INTEGER  I, J, K, ITERS, NUMTHR 

C  CLEAR  ALL  ARRAYS 

DO  100  I  -  I,  1000 

DN? (I)  “  0 
NPNCH(I)  -  0 
KENP(I)  =  0 
RRDHE(I)  -  0 
RRDDN(I)  **  0 
RRDDP(I)  -  0 
PFDHE(I)  -  0 
PFDDNC)  -  0 
FFDDP(I)  -  0 
PFTOT(I)  -  0 
P3REM ( I )  -  0 

i  T \  m  O 
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o  o  n  o  o  o 


FLOSS (I)  -  0 
DELTAP (I)  -  0 

100  CONTINUE 


C 

C  IMAXSQ  -  MAXIMUM  CURRENT  SQUARED  (AMPS) 

C  IMOPT  -  ASSUMED  VALUE  OF  MAXIMUM  ATTAINABLE  CURRENT  (AMPS) 

C  KTOPT  -  CORRESPONDING  PLASMA  TEMPERATURE  ASSUMING  T  GOES  AS  1**2  (KEV) 

C  CONST  -  CONSTANT  OF  PROPORTIONALITY  BETWEEN  TEMP  AND  CURRENT  (KEV/A**2) 

C  BOLTZ  -  BOLTZMANN'S  CONSTANT  (J/K) 

C  RHOI  -  INITIAL  FILL  GAS  DENSITY  (KG/M* *3) 

C  VRUN  -  PLASMA  SHEATH  RUNDOWN  VELOCITY  AT  THE  END  OF  THE  ANODE  (M/S) 

C  RA  ANODE  RADIUS  (M) 

C  RC  CATHODE  RADIUS  (M) 

C  LANODE  -  ANODE  LENGTH  (M) 

C  PNCHRAD  «  RADIUS  OF  PINCH  (M) 

C  LPNCE  -  PINCH  LENGTH  (M) 

C  FSNPLW  -  SNOWPLOW  EFFICIENCY  FACTOR,  FRACTION  OF  INITIAL  FILL  GAS 
C  WHICH  IS  ENTRAINED  IN  THE  RUNDOWN. 

C  FPNC3  -  PINCH  EFFICIENCY  FACTOR,  FRACTION  OF  GAS  IN  RUNDOWN  WHICH 
C  IS  TRAPPED  INSIDE  THE  PINCH. 

C  F  TOTALL  EFFICIENCY  -  FSNPLW*FPNCH 

C  DFRACT  -  PERCENTAGE  OF  DEUTERIUM  IN  FILL  GAS 

C  EEFRACT  -  PERCENTAGE  OF  HELIUM  IN  FILL  GAS 

C  AVMASS  -  AVERAGE  MASS  OF  PARTICLES  TRAPPED  IN  PINCE  (KG) 

C  ITERS  ”  NUMBER  OF  ITERATIONS  PERFORMED  DURING  EACH  PINCH 

C  NPINCH (I)  -  PINCH  NUMBER  DENSITY  FOR  Ith  ITERATION  (M**-3) 

C  REPRATE  -  NUMBER  OF  FIRINGS  PER  UNIT  TIME  (S**-l) 

C  PNCHTIM  -  DURATION  OF  PINCH  FORMATION  (S) 

C  KT  -  ENERGY  OF  PARTICLES  IN  PINCH  (KEV) 

I  IVOL  -  INITIAL  VOLUME  BETWEEN  ANODE  AND  CATHODE  (M**3) 

C  PNCHTIM  “  DURATION  OF  STABLE  PINCH  PHASE  (0) 

C  PVOL  -  FINAL  VOLUME  OF  PINCH  (M**3) 

C  PFDHE(I)  -  FUSION  POWER  FROM  D-HE3  REACTION  FOR  Ith  ITERATION  (W) 

C  PFDDN(I)  -  FUSION  POWER  FROM  DON  REACTION  FOR  Ith  ITERATION  (W) 

C  PFDDP(I)  -  FUSION  POWER  FROM  DDP  REACTION  FOR  Ith  ITERATION  (W) 

C  PFTOT(I)  -  TOTAL  FUSION  POWER  FOR  Ith  ITERATION  (PFDHE+PFDDN+PFDDP  (W) ] 

C  TPFDHE  -  TOTAL  FUSION  POWER  FROM  DHe3  REACTION  (W) 

C  TPFDDN  -  TOTAL  FUSION  POWER  FOR  DDn  REACTION  (W) 

C  TPFDDP  -  TOTAL  FUSION  POWER  FOR  DDp  REACTION  (W) 

C  TFP  -  TOTAL  FUSION  POWER  (W) 

C  PBREM(I)  -  RADIATIVE  LOSSES  DUE  TO  BREMSSTRAHLUNG  RADIATION  (W) 

C  PCYC(I)  -  RADIATIVE  LOSSES  DUE  TO  CYCLOTRON  RADIATION  (W) 

C  CYCREFL  -  FRACTION  OF  CYCLOTRON  RADIATION  RETAINED  BY  PLASMA 
C  PIN  -  POWER  NECESSARY  FOR  THE  OPERATION  OF  THE  FOCUS  (W) 

C  BPNCH  -  MAGNETIC  FIELD  IN  PINCH  (DETERMINES  PCYC)  (T) 

C  PLOSS(I)  -  TOTAL  POWER  LOST  OR  REQUIRED  TO  OPERATE  DEVICE  (W) 

C  DELTAP (I)  -  NET  POWER  INCREASE  OR  DECEREASE  (W) 

C  TPBREM  -  TOTAL  BREMSSTRAHLUNG  RADIATION  (W) 

C  TPCYC  -  TOTAL  CYCLOTRON  PADIATIC.J  GENERATED  (W) 

C  TPLOSS  -  TOTAL  RADIATIVE  POWER  LOSSES  (W) 

C  TDELTAP  -  TOTAL  NET  CHANGE  IN  POWER  (W) 

C  WCAP  -  INITIAL  ENERGY  STORED  IN  CAPACITOR  BANKS  (J) 

C  VOLT  -  CHARGING  POTENTIAL  OF  CAPACITOR  BANKS  (V) 

C  CAP  -  xNITIAL  EXTERNAL  CAPACITANCE  (CAPACITOR  BANK)  (F) 

C  LINIT  -  INITIAL  INDUCTANCE  OF  EXTERNAL  CIRCUIT  (H) 

C  LDOT  -  TIME  RATE  OF  CHANGE  OF  COAXIAL  INDUCTANCE  (H/S) 

I DOT  -  RATE  OF  CHANGE  OF  CURRENT  (A/S) 

A1-A4  “  CURVE  FIT  VALUES  TO  FIND  REACTION  RATE  PARAMETERS 
SIGVDHE  -  REACTION  RATE  PARAMETER  FOR  D-HE3  (M**3/S) 

StGVDQN  -  REACTION  RATE  PARAMETER  FOR  DDN  (M**3/S) 

SIGVDDP  -  REACTION  RATE  PARAMETER  FOR  DDP  (M**3/S) 

WDHE  -  ENERGY  RELEASED  PER  D-HE3  REACTION  (J) 

WDDN  -  ENERGY  RELEASED  PER  DDN  REACTION  (J)  89 
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C  WDDP  -  ENERGY  RELEASED  PER  DDP  REACTION  (J) 

C  DNP(I)  -  NUMBER  DENSITY  OF  DEUTERIUM  IN  PINCH  {K**~3) 

C  HENP  (I)  *»  NUMBER  DENSITY  OF  HELIUM  IN  PINCH  (M**~3) 

C  RRDHE(I)  -  RECTION  RATE  FOR  D-HE3  REACTION  <M**-3S**-1) 

C  RRDDN(I)  -  REACTION  RATE  FOR  DDN  REACTION  (M**-3S**-1) 

Z  RRDDP(I)  -  REACTION  RATE  FOR  DDP  REACTION  (M**-3S**-1) 

C  DSCBRG  -  TIME  FOR  FILL  GAS  TO  BE  DISCHARGED  (S) 

C  XSECTAR  -  CROSS  SECTIONAL  AREA  OF  FOCUS  DEVICE  (M**2) 

C  NFTHRST  -  THRUST  DUE  TO  EXPELLED  (NON-PINCH)  GASES  «*) 

C  FRACT  -  FRACTION  OF  ESCAPING  CYCLOTRON  RADIATION  ABSORBED 

C  IN  THE  WALLS  OF  THE  MIXING  CHAMBER  AND  ELECTRODES 

C  PA3SW  -  CYCLOTRON  POWER  ABSORBED  IN  THE  WALLS  6  ELECTRODES  ftfW) 

C  RHOCU  -  ELECTRICAL  RESISTIVITY  OF  COPPER  (OHM  M) 

C  PI2R  -  POWER  GENERATED  DUE  TO  OHMIC  HEATING  IN  THE  ELECTRODES  (MW) 

C  MH2  -  MASS  OF  DIATOMIC  HYDROGEN  MOLECULE  (KG) 

C  WASTE  -  WASTE  HEAT  DUE  TO  OHMIC  HEATING  AND  RADIATION  ABSORBED  (MW) 

C  MSFLW  -  COOLANT  MASS  FLOW  RATE  REQUIRED  TO  COOfrNKSTE  HEAT  AND 

C  KEEP  TURBINE  INLET  TEMPERATURE  LESS  THAN  2000  K  (KG/S) 

C  ADMSFLW  -  ANY  ADDITIONAL  HYDROGEN  USED  FOR  PROPELLANT  (KG/S) 

C  ENERGY  -  ELECTRICAL  ENERGY  PRODUCED  BY  TURBINE  AT  20%  EFFICIENCY  (MW) 
C  TAVG  -  MASS  AVERAGED  TEMPERATURE  OF  COOLANT  FROM  TURBINE  AND  ANY 

C  ADDITIONAL  PROPELLANT  FLOW  (K) 

C 

C  ISP  -  SPECIFIC  IMPULSE  (S) 

C  MUNOT  -  PERMITTIVITY  OF  FREE  SPACE  (H/M) 

C  GRAV  -  ACCELERATION  OF  GRAVITY  AT  EARTH'S  SURFACE  (M/'S* *2) 

C  TOTHRST  -  TOTAL  OF  ALL  THRUSTS  (N) , (L3F) 


★  A  ■*  + ■*■*++++  +  +*  +  +  + irir  it*  irititlritlrlr+’itit  frit +++'*+ +  * 

*******************  DEFINE  VALUES  OF  ALL  CONSTANTS  #****»»**»***»■»*•*** 

*  ********************************************************  ****■#•*  *»*•*»*#** 


DATA  MUNOT, PI, GRAV/ 1 . 257E-6, 3 . 1415, 9 . 8/ 

DATA  CYCREFL, DSCHRG/0 .6,1. OE-7/ 

DATA  DHEl,DHE2,DHE3,DHE4/0. 35715, -3. 32451, 10. 11363, -25. 66533/ 
DATA  DDNl,DDN2,DDN3,DDN4/0. 29811, -2. 08296, 5. 70135, -22. 0878/ 
DATA  DDPl,DDP2,DDP3,DDP4/0. 30795, -2. 12009, 5. 63718, -22. 03746/ 
DATA  WDHE, WDDN, WDDP, IMAGNET/2 . 93E-12, 5 . 24E-13, 6 .46E-13, 3 . 18E5/ 
DATA  MH2, BOLT Z, RHOCU, FRACT/ 3 . 34E-27, 1 .38E-23, 1 . 673E-8, 0 .2/ 

DATA  CPH2, CPELEC, CPHION, TDISION/4157, 1 . 517E7, 8267,5000/ 


Q  ******************************************************************** 

q  ***********  OPEN  INPUT  AND  OUTPUT  FILES  ************************** 

Q  *******♦*■>, ********************************************************** 

OPEN  (UNIT-1,  FILE-' test.dat' , STATUS-' OLD' ) 

OPEN  (UNIT-2,  FILE-' test. out' , STATUS-' NEW' ) 

OPEN  (UNIT-13, FILE-' twisp. out ', STATUS-' NEW' ) 

OPEN  (UNIT-14, FILE-' tw. out' , STATUS-' NEW' ) 

C  OPEN  (UNIT-15,  FILE-' isp. out' , STATUS-' NEW' ) 

READ(1,*)  90 

READ(1,*)  RA 

READd,*) 


READ ( I ,  *  )  RC 

READd,*) 

READ(1,*)  LANODE 
READd,*) 

READ(1,*)  REOI 
READd,*) 

R£AD(1,*)  VOLT 
READd,*) 

READd,*)  CAP 
READd,*) 

READd,*)  LINIT 
R£AD(1,*) 

READ(1,*)  SPECEN 
READ{1,*) 

READ (1, *)  FSKPLW 
READ(1/*) 

R£AD(1,*>  FPNCH 
READ(1,*) 

READ  (1,*)  LPNCH 
R£AD(1,*) 
R£AD(1,*)  PNCHRAD 
READ(1,*) 

R£AD(1, *)  DFRACT 
READd,*) 

READ (I#  *)  HEFRACT 
R£AD(1,*) 
READd,*)  REPRATE 
READ (1,  *) 
READ(1,*)  PNCBTIM 
READ(1,*) 
READ(1,*)  ITERS 
READ{1,*) 
READ(1,*)  DELTAV 
READ(1,*) 
READd,*)  NUMTHR 

CLOSE  (UNIT-1) 


Q  ************************************************************ 

C  *************  CALCULATE  FOCUS  PARAMETERS  ******************* 

g  *************************** ************************ ********* 

c 

MIN-1. 0E7 

DO  77  IMOPT-MIN, 1 ,5E7, 5 .0E6  " 

I  K-0 

DO  88  ADMSFLW-0 .0,30.0,0.5 

i  . 

XSECTAR  -  PI* (RC**2  -  RA**2) 

IVOL  -  XS ECT AR* LANODE 
PVOL  -  PI*PNCHRAD**2*LPNCH 
'i  F  -  FSNPLW  *  FPNCH 

AVMASS  -  HEFRACT* 5 . 0E-27  +  DFRACT*3 . 34E-27 
C  *******  MAXIMUM  CURRENT  ******** 

C  GIVEN  BY  DOLAN, "VOL  2 

IMAXSQ  -  2 . 704*SQRT ( (CAP*VOLT**3) / (MUNOT*LINIT*LOG (RC/RA) ) ) * 

£  ( (RA**2*RHOI/  MUNOT) **0 . 25)  . 
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o  o  ooo  o  n  oooo  Ov.oo 


IMAX  -  SQRT(IMAXSQ) 


*************  plasma  temperature  ********************* 


ASSUME  TEAT  PLASMA  TEMPE.  “T,E  " GALES  .*S  I  SQUARED 
ASSUME  MAXIMUM  ATTAINABLE  C.  -  TT  IS  IMOPT 

KT  -  MUNOT*IMAXSQ*AVMASS*LPNC3/(78.96*F*REOI*LANODE* 
£  (RC**2  -  RA**2) *1 . 6E-16) 

CONST  -  XT/IMAXSQ 

KTOPT  -  C0NST*IM0PT**2 


*******  RUNDOWN  VELOCITY  ******** 

TEE  RUNDOWN  VELOCITY  IS  CALCULATED  USING  TEE  MOMENTUM  EQUATION 
VRUN  -  SQRT (MUNOT * IMOPT * * 2 /  (78 . 96*RA**2*RHOI) ) 


*******  INITIAL  PINCH  NUMBER  DENSITIES  *•»***♦★** 

NPNCH (1)  -  F*RHOI*LANCDE*  (RC**2-RA**2) / (AVMASS*?NCERAD**2 
&  *L?NCH) 

DNP(l)  -  DFRACT*NPNCS(1) 

HENP(l)  -  EEFRACT *NPNC3 ( 1 ) 


**************************************************************  .-*•* 
***************  CALCULATE  REACTION  RATES  ************************ 

************************************************************** *y* 

X  -  LOGIO (KTOPT) 

SI6VDHE  -  1.0E-6M10**(DHE1*(X**3)+DEE2*(X*»2) 

&  tDHE3*X+DKE4)  )*1.5 

SIGVDDN  -  1 . OE-6* (10** (DDN1* (X**3) +DDN2* (X^*2) 

£  +DDN3*X+DuN4) ) 

SIGVDDP  -  1 .OE-6* (10** (DDP1* (X**3)+DDP2* (X**2) 
fi  +DDP3*X+DDP4) ) 

DO  99  1-1, ITERS 

RRDBE(I)  -  DNP (I) * (HEN? (I) *8X67085) 

RRDDN(I)  -  0. 5*DNP (I) * (DNP (I) *SIGVDDN) 

RRDDP(I)  -  0 . 5*DNP (I) * (DNP (I) *SIGVDDP) 


***********  DETERMINE  CHARGED  FUSION  POWER  FROM  PINCH  ****************** 

PFDHE(I)-RRDHEd) *WDKE*?VOL*REPRATE*PNCBTIM* (1.0/ITERS) *1.0E-6 
PFDDP  (I)-RRDDP  (I)  *WDDP*PVOL*RZPRATE*PNCHTIM*  (1 . 0/ITSRS)  *1 .OE-6 
PFDDN (I) -0 . 25*RRDDN (I) *WDDN*?VOL*R£?RATS*?NCHTIM*l . OZ-6 
fi  Ml.  0/ITERS) 

PFTOT (I)-PFDKE (I) +PFDDN (I) +PFDD? (I) 
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non 


***********  )  rTFPMTW:  HET  POWER  CHANGE  ***************************** 


VMAX-VCLT* ( (IMOPT+IMAGNET) /TMAX) 

PIN  -  VMAX  *  IMOPT  *  REPRATE  *  DSCHRG  *  1.0E-6 

C  BPNCH  IS  THE  MAGNETIC  FIELD  AT  THE  SURFACE  OF  THE  PINCH 
C  WHICH  IS  RESPONSIBLE  FOR  THE  CYCLOTRON  RADIATION  EMITTED 

BPNCH  -  MUNOT*F*IMOPT/ (2*PI*PNCHRAD) 

C  POWER  LOST  IS  THE  STEADY  STATE  LOSS  (PER  UNIT  VOLUME)  TIMES  THE 
C  PLASMA  VOLUME  TIMES  THE  TIME  PER  PINCH  TIMES  THE  NUMBER  OF  PINCHES 
C  PER  SECOND. 


PCYC (I) -6 .21E-17/ (8*PI) *MUNOT**2* (F*IMOPT) **2*LPNCH*NPNCH (I) * 

&  KTOPT* (l+KTOPT/146) *R£PRATE*PNCBTIM* (I . 0/ITERS) *1.0E-6 

PBREM(I)-5.35E-37*(F**2) *NPNCH (I) * (DNP (I)  +  (4*HENP  (I)  )  ) 

&  *  (SQRT  (KTOPT)  )  *PVOL*REPRATE*PNCHTIM*  (1 . 0/ ITERS) 

&  *1.0E-06 

ASSUME  THAT  CYCREFL  %  OF  CYCLOTRON  RADIATION  IS  RETAINED  IN  PLASMA 
PLOSS(I)  -  PBREM(I)  +  (1-CYCREFL) *PCYC (I) 

DELTAP (I)  -  PFTOT(I)  -  PLOSS(I) 


C 

C  ***********  DETERMINE  THE  THRUST  FROM  EXPELLED  FILL  GASES  ************ 
c  *******  ASSUMING  GASES  ARE  EXPELLED  FROM  THE  DEVICE  IN  IE-7  S  ******** 

NFTHRST  -  RHOI*VR.UN**2*XSECTAR*REPRATE*DSCHRG 


DNP  (I+U-DNP  (I)  -(RRDHE(I)  +2*RRDDN  (I)  +2*RRDDP  (I)  )  *PNCHTIM 
&  * (1.0/ITERS) 

HENP (1+1) -HFNP (I) -RRDHE(I) *PNCHTIM* (1 . 0/ITERS) 

NPNCH (1+1)  -  DNP ( 1+1)  +  HENP (T+l) 

99  CONTINUE 

TFP  -  0 
TPLOSS  -  0 
TDELTAP  -  0 
TPFDHE  -  0 
TPFDDN  -  0 
TPFDDP  -  0 
TPBREM  -  0 
TPCYC  -  0 


C  TOTAL  ELEMENTS  IN  ALL  ARRAYS 
C 

DO  200  J  -  1, ITERS  93 


TFP  -  TF?  +  PFTOT(J) 


o  o  o  o  o  no  o  o  no 


TP LOSS  -  TP LOSS  +  PLOSS(J) 

TDELTAP  -  TDELTAP  +  DELTAP (J) 

TPFDHE  -  TPFDHE  +  PFDHE(J) 

TPFDDN  -  TPFDDN  +  PFODN(J) 

TPFDDP  -  TPFDDP  +  PFDDP(J) 

TPBREM  -  TPBREM  +  PBRZM(J) 

TP CYC  -  TPCYC  +  PCYC(J) 

200  CONTINUE 

C  *******  CYCLOTRON  RADIATION  ABSORBED  IN  WALL  AND  ELECTRODES  ******* 
PABSW  -  FRACT  *  (1-CYCREFL) *TPCYC 

C  *******  TOTAL  POWER  DISSIPATED  IN  ELECTRODES  BY  OHMIC  HEATING  ******* 

PI2R  -  mOPT**2*RHOCU*LANODE*(l/(PI*RA**2)+l/(RI*(RC**2 
&  -RA**2 ) ) ) *1 . 0E-6*DSCHRG*REPRATE 

c  *************  POWER  TO  BE  REMOVED  FROM  MAGNET  (MW)  **********>***** 

P MAGNET  -  0.01 

c  *******  TOTAL  POWER  TO  BE  REMOVED  FROM  THE  WALLS  AND  ELECTRODES  **** 

C  ***********************  ajid  MAGNET  ************+******■**'***’*■*'**'**'*—* 

WASTE  -  PI2R  +  PABSW  +  PMAGNET 


**  NEED  TO  KEEP  INLET  TEMPERATURE  TO  TURBINE  BELOW  ABOUT  2000K  ***** 
*****MASS  FLOW  REQUIRED  TO  DO  SO  IS  GIVEN  BY  ****** 

MSFLW  -  WASTE*61.14/(2000  -  20) 

********  ASSUME  TURBINE  TO  BE  20%  EFFICIENT  ************ 

ELECTRICAL  POWER  GENERATED  IS  20%  OF  WASTE  HEAT 

ENERGY  -  0 . 2*WASTE 

**»  ASSUME  TURBINE  EXIT  TEMPERATURE  IS  ABOUT  700K,  WHERE  IT  "CAN 
BE  MIXED  WITH  ADDITIONAL  MASS  FLOW  AT  20K  ***** 


TAVG  -  (MSF LW*700  +  ADMSFLW*20) / (ADMSFLW+MSFLW) 


****  ASSUME  GAS  ABSORBS  MOST  FUSION  POWER  PRODUCED,  MIXES  UNIFORMLY- 
AND  EXITS  AT  A  UNIFORM  TEMPERATURE  OF  TOUT  ****** 

ASSUME  GAS  ABSORBS  HEAT  UP  TO  5000K  AS  H2  WITH  GAMMA-1.40 
THE  HEAT  REMOVED  FROM  THE  SYSTEM  IN  DOING  SO  IS: 


Q  ***** 

QDOTREM- (ADMSFLW+MSFLW) *CPH2* (TDISION-TAVG) 

£  ***** 

TEMP-TAVG+1 .0E6*TDELTAP/ (CPH2* (ADMSFLW+MSFLW) ) 
IF  (TEMP  .LT.  TDISION)  THEN 
GOTO  83 
END  IF 


C  THIS  LEAVES  A  TOTAL  OF  QLEFT  TO  BE  ABSORBED  BY  A  COMPLETELY 
C  DISSOCIATED  AND  IONIZED  COMBINATION  OF  AN  ELECTRON  GAS 
C  AND  A  GAS  OF  HYDROGEN  IONS.  94 


OOOO  O  O  OO  O  O  OO  OoOO 


^  TC  X  *  X 

Q LEFT-1 . 0E6  *TDELTAP-QDOTRZM 

C  ’*»*** 

C  NOW  THE  POWER  IS  ABSORBED  BY  AN  ELECTRON  GAS  AND  A  HYDROGEN 
ION  GAS.  TEE  CONSTANT  VOLUME  HEAT  CAPACITY  OF  A  FREE 
C  MONATOMIC  GAS  CAN  BE  FOUND  USING  FERMI -DIRAC  STATISTICS  TO  BE 
.  C  CV-(N/2)*R,  WHERE  R  IS  TEE  UNIVERSAL  GAS  CONSTANT  AND  N  IS  THE 

C  NUMBER  OF  DEGREES  OF  FREEDOM  OF  EACH  PARTICLE,  IN  OUR  CASE,  3. 

Q  ****************  H  aa  8.3143  J/incl  X  ******************* 

C  ASSUMING  AN  IDEAL  GAS,  CP  -  CV  +  R 
C  ie.  C?  -  (5/2) *R  -  20.786  J/mol  X 
C  1  MOLE  OF  ELECTRONS  IS  5.48E-7  XG 
C  1  MOLE  OF  H+  IONS  IS  1.006E-3  XG 
C  CP ELEC  -  1.517E7  J/XG  X 

C  CPHION  -  8267  J/XG  X 
C 

C  QLEFT  -  MOOT*  (CPELEC+CPBION)  *DELTAT 

C  AFTER  ABSORPTION  OF  THIS  ENERGY,  THE  PLASMA  IS  ASSUMED  TO 
C  COME  TO  THERMAL  EQUILIBRIUM  AT  STAGNATION  CONDITIONS,  ie.  V-0. 


Q  ***** 

T ST AGE-TD I S ION  +  QLEFT/ ( (ADMSFLW+MSFLW)  *  (CPHION*0 . 999455) ) 
TSTAGE-TDISION  +  QLEFT/  ( (ADMSFLW+MSFLW)  *  (C?£L£C*0 .000545) ) 

Q  ***** 

ASSUME  TEAT  THE  FLOW  NCW  ENTERS  A  MERIDIONAL  MAGNETIC  NOZZLE 
WITH  ONE  COIL  AT  THE  THROAT  OF  TEE  NOZZLE. 

THE  SPECS  FOR  SUCH  A  NOZZLE  ARE  GIVEN  IN  THE  AL  REPORT 
"CHARACTERIZATION  OF  .PLASMA  FLOW  THROUGH  MAGNETIC  NOZZLES". 

ACCORDING  TO  TEE  SPECS,  TSTAG/TTHR  -  1.35 
AND  VEXIT/VTHROAT  -  2.0 

***** 

TTHRTH  -  TSTAGH/1 .35 
TTHRTE  -  TSTAGE/1 .35 

***** 


FROM  CONSERVATION  OF  ENERGY:  CP*DELTAT- (1/2)  *VTHROAT**2 
THERMAL  ENERGY  IS  CONVERTED  TO  ENTHALPY  OF  THE  PLASMA 

***** 

VHIONTH-SQRT (2*CPHION* (TSTAGE-TTHRTH) ) 

VELECTH-SQRT (2*CPELEC* (TSTAGE-TTERTE) ) 

***** 


FLOW  EXITS  TWICE  AS  FAST  BECAUSE  OF  EXPANSION 
THROUGH  THE- NOZZLE.  ' .  ...... 

VHIONEX  -  2 . Q*VHICNTH 
VELECEX  -  2 . 0*VELECTH 


C  RESULTING  THRUST  FROM  PLASMA 

* 

£  ***** 

ELECTER  -  5 . 45E-4* (ADMSFLW+MSFLW) "VELECEX 
HIONTHR  -  0.999455* (ADMSFLW+MSFLW) "VHIONEX 

95. 


boooowo  oooo  n  o  ooooooo  o  o  o  o  o  ooooo  o  o  oono 


PROPTER  -  ELECTSR  +  HICNTER 
TOTERST  -  (NFTSRST  +  PRCPTER) 

IS? -TOT EES 7/  (GRAV*  (RSOI*IVOL*RZ?RATE*  { 1-? )  +ADMSFLW+MSFLW) ) 
VZX-IS?*GRAV 

Q 

*******  ESTIMATE  masses  for  mars  mission  ******************** 

ASSUME  FIXED  PAYLOAD  MASS  OF  ABOUT  100  METRIC  TONS 

*  *  if  ★★ 

MPAYLD  -  1.0E5 

*  *★  it* 

DETERMINE  TEE  NECESSARY  CAPACITOR  MASS  USING  INPUTTED 

SPECIFIC  ENERC-Y 

ELECTRIC  ENERGY  FROM  CAPACITOR  BANKS 


***** 

WNOT-O . 5*CAP*VOLT**2 

ELECTEN  -  WNOT*  { (IMOPT+IMAGNET)  /XMAX)  **2 . 0 
***** 


NEED  ELECTEN  IN  JcJ  AND  SPECEN  IN  kJ/kg 


***** 

MCA?  -  NUMTER*ELECTEN/1*'00/S?ECEN 
***** 


MASS  OF  MAGNET  TO  BE  USED  AT  TEE  CENTER  OF  TEE  MAGNETIC  NOZZLE 
ASSUME  A  FIELD  OF  2  TESLA  IS  NEEDED  AT  TEE  TEROAT  WHICE  EAS 
A  RADIUS  OF  10  CM.  THE  COPPER  MAGNET  WILL  THUS  EAVE  AN  INNER  RADIO'S 
CF  10  CM  AND  AN  OUTER  RADIUS  CHOSEN  AS  50  CM  IN  ORDER  TO 
MINIMIZE  RESISTANCE.  THE  MAGNET  WILL  BE  PULSED  100  TIMES  PER 
SECOND  AND  EACH  PULSE  WILL  LAST  A30UT  10**-4  SECONDS.  TEE  LENGTH 
OF  TEE  MAGNET  WAS  CHOSEN  TO  BE  1  CM. 

***** 

MMAGNST  -  67 . 55*NUMTER 
***** 


POWER  DISSIPATED  IN  TEE  MAGNET  IS  I**2*R*?NCSTIM*RE?RAT£ 

OR  ABOUT  10  KW.  NOTE  TEAT  TSESE  NUMBERS  GIVE  ENERGY  DENSITIES 
IN  TEE  MAGNET  WSXCH  ARE  MUCH  LESS  TEAN  TEE  MAXIMUM  TOLERABLE 
ENERGY  DENSITY  AT  WHICH  COPPER  BEGINS  TO  MELT. 


******  CALCULATE  TEE  AMOUNT  OF  PROPELLANT  NEEDED  FROM  THE 
MISSION  DELTA  V  AND  TEE  EXHAUST  VELOCITY  USING  TEE 
ROCKET  EQUATION. 

AS  A  FIRST  APPROX.  TO  THE  INITIAL  MASS,  CONSIDER  ONLY 
PAYLOAD,  CAPACITORS,  AND  MAGNET  MASSES.  THIS  MAKES  TEE 
CALCULATION  CONSIDERABLY  EASIER  AND  INTRODUCES  ONLY  A 


on  o  o  non  o  o  non  noono 


c 

MPROP-  (EXP  (DELTAV/VEX)  -1)  *  (MPAYLD+MCAP+MMAGNET) 
MPROSTR-O . 15*MPR0P 

•  ******  total  born  time  ************ 

TBURN  -  MPROP/ (ADMSFLW+MSFLW) 

C  MASS  OF  DEUTERIUM  AND  HELIUM  3  FUEL 
«  £  ***** 

MFUEL  -  NUMTHR*RHOI*IVOL*R£PRATE*TBURN 
MFUELSY  -  0 . 1 *MFUEL 


CALCULATE  MASS  OF  SHIELD  NECESSARY  TO  KEEP  NEUTRON  FLUENCE 
BELOW  10**13  FOR  A  MISSION  THAT  IS  DAYS  LONG 

LITHIUM  HYDRIDE  SHIELD  THICKNESS  IS  IN  METERS 


ASSUME  SHIELD  IS  1  METER  FROM  NEUTRON  SOURCE  AND  SUBTENDS 
AN  ANGLE  SUCH  THAT  ABOUT  12.5%  OF  ALL  NEUTRONS  RELEASED 
IN  THE  DDn  REACTION  HIT  THE  SHIELD 

★  *  ★  ★  * 

THICK  -  0  .  l*LOG(Q  .  125*RRDDN (ITERS+1)  *  PNCHTIM*  REP  RATE  *P  VOL 
&  *T8URN/ (8 6400* ( 4 *P 1*1. 157E12) ) ) 


ASSUME  SHIELD  HAS  A  CROSS  SECTIONAL  AREA  OF  ONE  METER 
MASS  -  DENS IT Y* AREA*THICKNESS  (DENSITY  OF  LiH  IS 
APPROXIMATELY  725  KG/M**3) 


*  *  it  -k  it 

MSHIELD  -  NUMTHR*725.0*PI*THICK 

★  ir  tt  *  if 


****************  TOTAL  MASS  CALCULATION  (IN  KG)  *************** 


MTOT  -  MPAYLD+MCAP+MPROP+MPROSTR+MFUEL+MFUELSY 
&  +MSHIELD+MMAGNET 

WRITE (13, *)  ISP , NUMTHR*TOTHRST/ (GRAV*MTOT) 

WRITE  (14,  *)  MSFLW+ADMSFLW, NUMTHR*TOTHRST/  (GRAV*MTOT) 
C  WRITE (15,*)  MSFLW+ADMSFLW, ISP 
K-K+l 

88  CONTINUE 
PRINT*, K 
77  CONTINUE 


C  **************  WRITE  RESULTS  TO  OUTPUT  FILE  IMP. OUT  ****************** 

999  WRITE (2,*)  'INITIAL  FILL  GAS  DENSITY' 

WRITE (2,*)  RHOI 

WRITE (2,*)  'MAXIMUM  CURRENT  IN  AMPS' 

WRITE  (2,*)  IMAX  97  ■ 


WRITE (2,*) 
WRITE (2/ *) 
WRITE (2,  *) 
WRITE (2,  *) 
WRITE (2,*) 
WRITE (2,*) 
WRITE (2,  *) 
WRITE (2,  *) 
WRITE (2, *) 
WRITE (2,*) 
WRITE (2,  *) 
WRITE (2,*) 
WRITE (2, *) 
WRITE (2,*) 
WRITE (2, *) 
WRITE (2,*) 
WRITE (2,*) 
WRITE (2,  *) 
WRITE (2,*) 
WRITE (2,*) 
WRITE (2, *) 
WRITE (2,*) 
WRITE (2,*) 
WRITE (2,*) 


'PLASMA  PINCH  TEMPERATURE  IN  KEV' 

XT 

'AT  A  CURRENT  OF  (MA) ' 

IMOPT 

'PLASMA  PINCH  TEMPERATURE  IN  KEV' 

KTOPT 

'RUNDOWN  VELOCITY  AT  THE  END  OF  THE  ANODE:. IN  H/S' 
VRUN 

'INITIAL  AND  FINAL  D  NUMBER  DENSITY  IN  M**^3' 
DNP(l) ,DNP( ITERS +1) 

'INITIAL  AND  FINAL  HE  NUMBER  DENSITY  IN  M**-3' 

HENP ( 1 )  , HENP ( ITERS+I ) 

'FRACTION  OF  DEUTERIUM  BURNED' 

1- (DNP (ITERS+1) /DNP (I) ) 

'FRACTION  OF  HELIUM-3  BURNED' 

I- (HENP (ITERS+1) /HENP (1) ) 

'REACTION  RATE  PARAMETERS  FOR  DHe3,  DDn,  AND  DDp  ' 
'IN  M**3/S' 

SIGVDHE 

SIGVDDN 

SIGVDDP 


WRITE (2,*)  'INITIAL  AND  FINAL  REACTION  RATES  FOR  DH«3,  DDn, 
&  AND  DDp  IN  M»*-3S**-1' 


WRITE (2,*)  RRDHE(l) ,RRDHE (ITERS) 

WRITE (2,  *)  RRDDN(l) ,RRDDN (ITERS) 

WRITE (2, *)  RRDDP(l) ,RRDDP( ITERS) 

WRITE (2,*) 

WRITE (2,  *)  'CHARGED  PARTICLE  FUSION  POWER  FROM  DH«3,  DDn,  AND 
fiDDp  IN  MEGAWATTS' 


WRITE (2,*) 
WRITE (2,*) 
WRITE (2 ,  *) 
WRITE (2,*) 
WRITE (2,*) 
WRITE (2,*) 
WRITE (2,*) 
WRITE (2,*) 
WRITE (2,*) 
WRITE (2,*) 


TPFDHE,TPFDDN, TPFDDP 

'TOTAL  FUSION  POWER  IN  MEGAWATTS' 

TFP 

'POWER  NEEDED  TO  OPERATE  FOCUS  IN  MEGAWATTS' 

PIN 

'BREMSSTRAHLUNG  AND  CYCLOTRON  LOSSES  IN  MEGAWATTS' 

TP3REM, (1-CYCREFL) *TPCYC 

'TOTAL  POWER  LOSSES  IN  MEGAWATTS' 

TPLOSS 


IF  (TDELTAP-PIN  .LT.  0)  THEN 

WRITE (2,*)  'NET  DECREASE  IN  POWER  IN  MEGAWATTS' 

WRITE (2,*)  TDELTAP-PIN 
ELSEIF  (TDELTAP-PIN  .GT.  0)  THEN 

WRITE (2,*)  'NET  INCREASE  IN  POWER  IN  MEGAWATTS' 

WRITE (2, *)  TDELT AP -P IN -PMAGNET 
ELSE 

WRITE (2,*)  'THERE  IS  NO  NET  CHANGE  IN  REACTOR  POWER' 
END  IF 


WRITE (2,*) 
WRITE (2,*) 
WRITE (2,*) 
WRITE (2,*) 
WRITE (2,*) 
WRITE (2,*) 
WRITE (2,*) 
WRITE (2, *) 
WRITE (2,*) 


'TOTAL  ELECTRICAL  POWER  PRODUCED  IN  MEGAWATTS' 
ENERGY 

'TOTAL  MASS  FLOW  RATE  IN  KG/S' 

ADMSFLW+MSFLW 

'ION  STAGNATION  TEMPERATURE  IN  K,eV' 

TSTAGH, TSTAGH/12000 

'ELECTRON  STAGNATION  TEMPERATURE  IN  K,eV' 
TSTAGE , TSTAGE/ 12000 

'tow  1NH  ft rpTQON  FVTT  VELOCITIES  IN  M/S' 


WRITE (2,*)  VHIONEX, VELECEX 

WRITE  (2,  *)  'FINAL  PROPELLANT  THRUST  IN  N' 

WRITE (2,*)  PROPTHR 

WRITE (2,*)  'TOTAL  BURN  TIME  IN  S, DAIS' 

WRITE (2,*)  TBURN,TBURN/86400 

WRITE (2,  *) 

WRITE (2,*)  *************  SYSTEM  MASSES  IN  KG' 

WRITE (2,*)  'PAYLOAD  MASS : ' ,MPAYLD 
WRITE (2,*)  'PROPELLANT  MASS:',MPROP 

WRITE (2, *)  ' PROPELLANT  SYSTEM  AND  STRUCTURE : ' ,MPROSTR 

WRITE (2,  *)  'FUEL  MASS:',MFUEL 

WRITE (2,*)  'FUEL  SYSTEM  MASS : ' , MFUELSY 

WRITE (2,*)  'ELECTRIC  ENERGY  FROM  CAPACITORS  (J) ' , ELECTEN 

WRITE ( 2 , * )  '  CAPACITOR  MASS : ' , MCAP 

WRITE (2,*)  'SHIELD  MASS : ' ,MSHIELD 

WRITE (2,*)  'MAGNET  MASS : ' ,MMAGNET 

WRITE (2,*)  ***********  TOTAL  MASS:', MTOT 

WRITE (2,  *) 

WRITE (2,*) 'TOTAL  THRUST  IN  NEWTONS, LBF' 

WRITE  (2,  *)  NUMTHR*TOTHRST,NUMTHR*TOTHRST/4 .4482 
WRITE (2,*) 

WRITE (2,*)  'THRUST  TO  WEIGHT  RATIO' 

WRITE  (2,  *)  NUMTHR*TOTHRST/  (MTOT*GRAV) 

WRITE (2,*)  'SPECIFIC  IMPULSE  IN  SECONDS' 

WRITE (2, *)  ISP 

WRITE (2, *)  ' ******************************************* 

WRITE (2,*) 

CLOSE (UNIT-2) 

CLOSE (UNIT-13 ) 

CLOSE (UNIT-14) 

CLOSE (UNIT-15) 

CLOSE (UNIT-16) 

CLOSE (UNIT-17) 

END 


99 
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